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Abstract 

In this paper, we consider the hamihonian formulation of nonholonomic systems with symme- 
tries and study several aspects of the geometry of their reduced almost Poisson brackets, including 
the integrability of their eharaeteristie distributions. Our starting point is establishing global for- 
mulas for the nonholonomic Jacobiators, before and after reduction, which are used to clarify the 
relationship between reduced nonholonomic brackets and twisted Poisson structures. For certain 
types of symmetries (generalizing the Chaplygin case), we obtain genuine Poisson structures on the 
reduced spaces and analyze situations in which the reduced nonholonomic brackets arise as gauge 
transformations of these Poisson structures. We illustrate our results in mechanical examples, and 
also show how to recover several well-known facts in the special ease of Chaplygin symmetries. 
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1 Introduction 



Nonholonomic systems are mechanical systems with constrains in their velocities that are not derived 
from constraints in the positions. Examples include rolling contact constraints and special types of 
sliding constraints, see e.g. [5} UH [6]. In this paper we use the hamiltonian viewpoint to nonholonomic 
systems to study several aspects of their geometry in the presence of symmetries. 

Nonholonomic systems are described by a configuration manifold Q, a lagrangian C : TQ — >■ M, 
assumed to be of mechanical type, and a non-integrable distribution D C TQ, defining the permitted 
velocities of the system. These systems admit a hamiltonian formulation (developed e.g. in [H [191 [25]) 
in terms of triples (7W, vr^h, ^ai); where A4 is the submanifold of T*Q given by the image of the 
constraint distribution D under the Legendre transform Leg : TQ —^T*Q, Hj^ is the restriction to 
M of the hamiltonian function on T*Q corresponding to C, and TTnh is a bivector field on Ai naturally 
defined from D and the canonical symplectic form Qq on T*Q. The associated nonholonomic bracket 
(see [ia[2a[29]) {■, -U is given by 

{f,9U = 7r^M,dg), f,geC^iM), 

and the evolution of the nonholonomic system is given by the flow of the nonholonomic vector field 

An important feature that distinguishes nonholonomic systems from usual hamiltonian systems is that 
the bracket {•, -jnh does not satisfy the Jacobi identity, so it is just an almost Poisson bracket. In 
fact, it can be shown that the characteristic distributioio of {•, -Inh is not integrable, in contrast with 

^the distribution on M generated by vector fields of the form Xf :— {■, /}nii, for / G C°°{A4). 
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Poisson structures, whose characteristic distributions are always tangent to symplectic leaves. Since 
almost Poisson brackets and bivector fields are equivalent, we will treat them indistinguishably. 

The present paper is concerned with the lack of integrability of nonholonomic brackets, especially 
in the case of nonholonomic systems with symmetries, i.e., when Q is equipped with a free and proper 
action of a Lie group G preserving D and C. In this case, M. inherits a G-action such that TTn,,, T-Lj^ 
and Xjih are G-invariant. As a consequence the orbit space M/G \s endowed with a bivector field 
7r"^'j, (defining the reduced nonholonomic bracket {■,-}'^^^, a reduced hamiltonian H^ed) and a reduced 
nonholonomic vector field defining the reduced dynamics and given by 



In many concrete situations, the reduced bivector tt'^^^ has better integrability properties than TTnh- 
Indeed, in some cases the reduced bracket is a Poisson bracket, possibly after a time reparametrization, 
see e.g. [Dj [HI [151 [H [El [13 [2Ql [27]. This is a useful scenario for the analysis of the integrability of 
the system [T5] or to develop a Hamilton- Jacobi theory |271 [T2] . 

Even if vrj* is not a Poisson structure, it may still have an integrable characteristic distribution, 
which means that it is described by almost symplectic leaves. Note that the existence of this foliation 
may give information about the nonholonomic dynamics, for example concerning the presence of con- 
served quantities, since each leaf is invariant by the flow of ^"Jj,. One of our motivations in this paper 
is understanding geometric features of reduced nonholonomic brackets, particularly the integrability of 
their characteristic distributions. 

It is worth noticing that one may have many brackets generating the nonholonomic vector field. 
More precisely, we will say that an almost Poisson bracket {•, •} on M (respectively on M./G) describes 
the dynamics if = X^i, (respectively {•,'Hred} = -^r"d)- So, even if {^-j^^d does not satisfy 

properties of interest, one can still hope of finding another bracket describing the dynamics with the 
desired properties. A general mechanism to produce brackets describing the dynamics was introduced 
in [2]. It is based on the fact that one can use 2-forms to modify almost Poisson brackets via the 
so-called gauge transformations [28]. In this procedure, 2-forms B on M are used to "deform" the 
bivector field 7r„h keeping its characteristic distribution unchanged. The resulting modified bracket will 
be denoted by tts, and two brackets related by a gauge transformation are called gauge related. If B is 
annihilated by then ttb still describes the dynamics. An interesting point is that the reduction of 
gauge related bivector fields need not be gauge related on M/G. In fact, one may verify (see [2]) that 
vTnh and -Kb may lead to reduced brackets vr^^'Ji, vr^^^ with fundamentally different geometric properties. 
For example, irf^^ can have integrable characteristic distribution (or even be a Poisson bracket) while 
7r°'Ji does not. So we will be concerned in this paper not only with the geometry of vrj'^'jj but also with 
other brackets of the form vr^^^ describing the reduced dynamics. 

In order to describe our results, it is convenient to recall some fundamental facts regarding Chaplygin 
symmetries (see e.g. |2H [5]). A nonholonomic system with symmetries is called Chaplygin if the 
constraint distribution D complements the vertical space V C TQ with respect to the group action: 



In this case the reduced nonholonomic bivector vr^^'jj is nondegenerate |2T], so it is equivalent to a 
nondegenerate 2-form i}'^^^ on M/G. The nonholonomic character of the reduced system is measured 
by the non-closedness of ^2"^'^. The following results are well known. 



nh 
red ■ 



TQ = D®V. 



(1.1) 
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(a) The differential of JlJ^^ is given by [5] 

d^t^ = -d{J,iq,^^, (1.2) 

where (^7, /C)rcd is a 2-form on Ai/G explicitly obtained from the canonical momentum map for 
the lifted action on T*Q and the curvature of the constraint distribution (viewed as a principal 
connection). 

(6) There is an identification |21j 

M/G ^T*{Q/G). (1.3) 
Hence M/G has a natural symplectic structure q/q, usually not describing the dynamics, 
(c) The relation between and ^q/g is given by [22] 

^t. = ^QJG-{J.K,),^,. (1.4) 

Note that, in particular, this explains ()1.2p in (a). 

In this paper, we provide extensions of the results in (a), (6), (c) to more general situations. Instead 
of assuming condition (jl.ip , we will replace it by the weaker dimension assumption [6] 

TQ = D + V. (1.5) 

Note that many classical nonholonomic systems with symmetries satisfy (jl.Sp but not (jl.ip . see e.g. 
[51 [3]. Some explicit examples will be discussed in the end of the paper. 

Our starting point is establishing a global, coordinate-free formula for the Jacobiator of the non- 
holonomic bracket {•, -j^h on M (extending the local formula given in [251 sec. 2.5]), not making use of 
symmetries. Once symmetries are present (satisfying ()1.5p ). the geometric meaning of this Jacobiator 
formula becomes more clear, and leads to an expression for the Jacobiator of the reduced nonholonomic 
bracket {-j-jrcd M/G. In fact, all our formulas can be easily extended to more general brackets 
obtained from {•, -jnh and {•, ■}'^^^ via gauge transformations by 2-forms. 

One application of our nonholonomic Jacobiator formulas is clarifying the link between reduced 
nonholonomic bivectors tt^^JJj (or, more generally, irf^^) and twisted Poisson brackets. This issue was 
already raised in e.g. [231 Sec. 1.8]. Recall that a twisted Poisson structure [28] is a special type of 
almost Poisson bracket whose Jacobiator is controlled by a closed 3-form; an important feature of 
twisted Poisson structures is that they always have integrable characteristic distributions. Prom the 
formula describing the reduced Jacobiators, we find explicit conditions for the reduced brackets to be 
twisted Poisson as well as a mechanical interpretation of the 3-form with respect to which the reduced 
bracket is twisted. Moreover, in the special case of Chaplygin symmetries, the Jacobiator of {•, ■}'f^^ 
can be expressed in terms of dOs^"^^ (see ()2.6p and ()2.1ip below) , and our formula recovers ()1.2p in (a) . 

In order to generalize {h) and (c), we consider symmetries satisfying (|1.5p and fix an invariant 
vertical complement W C y of D, so that TQ = D ® W , satisfying an additional condition saying 
that W can be realized as the vertical space for a subgroup of the group G of symmetries. Once such 
W is chosen, we prove that one naturally obtains a Poisson structure A on the reduced space Ai/G. 
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We also show that the submanifold W° /G of T*Q/G, where W° C T*Q is the annihilator of W, has a 
natural Poisson structure Aq, and that there is an identification 

~ W°/G 

which is a Poisson diffeomorphism. In the Chaplygin case, the vertical complement W is necessarily 
equal to V, the quotient V° /G is canonically identified with T*{Q/G), and Aq coincides with the 
Poisson structure defined by canonical symplectic 2-form ^q/g, thus recovering the result in (b). 

In general, the Poisson structure A on Ai/G does not describe the nonholonomic dynamics, so one 
is led to compare the Poisson structure A with the reduced nonholonomic bivector vr"^'^,, as done in (c). 
Note that, in terms of bivector fields, the relation in (|1.4|) is equivalent to saying that the bivectors 
corresponding to Q'^^^ and r^g/c are gauge related by the 2-form (J', /C)red- Unlike the Chaplygin case, 
we observe that, in general, vr^^'jj is not necessarily gauge related to A. However we obtain sufficient 
conditions, fulfilled in many examples, guaranteeing that tt^^^ (or, more generally, irf^^) is obtained 
from the Poisson bivector A by a gauge transformation. We show that this happens e.g. for the 
vertical rolling disk, the snakeboard and the Chaplygin ball. Notice that proving that the reduced 
dynamics is described by a bivector gauge related to the Poisson structure A has several implications; 
for example, the nonholonomic flow is restricted to the symplectic leaves of A, so Casimirs of A are 
conserved quantities of the system. 

Outline of the paper. In Section [2] we review the basics of almost Poisson structures, gauge transforma- 
tions and nonholonomic systems that are useful to our paper. In Section [3] we prove the nonholonomic 
Jacobiator formula (see Theorem 13.21 and Corollary 13. 3p . Nonholonomic systems with symmetries 
satisfying the dimension assumption (jl.Sp are considered in Section SI We prove a formula for the 
Jacobiator of nonholonomic brackets in this case (Theorem 14. 5p as well as a formula for the Jacobiator 
of the reduced brackets (Corollary 14. 7p : here we also establish the connection with twisted Poisson 
structures (Corollarv I4.10p . We show how the Jacobiator formulas simplify once additional conditions 
are imposed on symmetries (see Theorem 14. 131 and Corollarv I4.14p . and illustrate our results in the ex- 
ample of a nonholonomic particle and Chaplygin systems. In Section [5] we define the Poisson structure 
A on Ai/G, based on a suitable choice of G-invariant vertical complement W oi D (see Prop. 15. 2p . 
We also define, in two equivalent ways, a Poisson structure Aq on W°/G (see Propositions 15.61 and 15. SP 
which is independent of the kinetic metric, and prove that it is Poisson diffeomorphic to {Ai/G,A) 
(Prop. I5.10p . In Section [6] we give a description of the symplectic leaves of the Poisson structure A on 
M/G (in terms of the leaves of the canonical Poisson structure on T*Q/G, see Theorem 16. ip . which 
turns out to be related to the nonholonomic momentum map (see Theorem 16.51 and Corollarv 16. 6p . 
We also analyze, in Theorem 16.81 when the reduced nonholonomic brackets can be written as gauge 
transformations of the Poisson structure A. In Section [71 we work out several mechanical examples 
illustrating our results. We have also included an Appendix collecting some facts about the reduction 
of presymplectic forms used in the paper. 

Acknowledgments: I thank FAPERJ (Brazil) and the CMC Network (projects MTM2010-12116-E and 
MTM2012-34478, Spain), especially J.C. Marrero, E. Padron and D. Martin de Diego, for supporting 
this project. I have also benefited from attending the conferences Applied Dynamics and Geometric 
Mechanics, held in Oberwolfach, and the Focus Program on Geometry, Mechanics and Dynamics, the 
Legacy of Jerry Marsden, held at the Fields Institute, so I thank the various organizers, particularly A. 
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Bloch and T. Ratiu, for their invitations. I am grateful to Alejandro Cabrera, Oscar Fernandez, Luis 
Garcia-Naranjo, David Iglesias-Ponte, Jair Koiller and Tom Mestdag for several interesting discussions 
and encouragement, and Henrique Bursztyn for his many comments on the manuscript. I particularly 
thank Oscar Fernandez for pointing out the snakeboard example. 

2 Almost Poisson structures and nonholonomic systems 

2.1 Almost Poisson brackets and gauge transformations. 

In this section we recall basic definitions related to almost Poisson manifolds. 

An almost Poisson bracket {•, •} on a manifold P is an M-bilinear bracket {•, •} : C°°{P) x C°°{P) — )■ 
C°°{P) that is skew-symmetric and satisfies the Leibniz condition: 



If the Jacobi identity is also satisfied then the bracket {•, •} is called a Poisson structure. 

Let P be a manifold equipped with an almost Poisson bracket {•,•}. The hamiltonian vector field 
Xf oi a function / € C°°{P) is the vector field on P defined by 



for all g G C°°{P). The characteristic distribution of {•, •} is the distribution on the manifold P whose 
fibers are spanned by the hamiltonian vector fields. In general the characteristic distribution is not 
integrable. If it is integrable, so that it is tangent to leaves (of possibly varying dimensions), then each 
leaf inherits a nondegenerate 2-form; i.e., an almost Poisson structure with integrable characteristic 
distribution gives rise to a (singular) foliation with almost symplectic leaves. When an almost Poisson 
bracket is Poisson then its characteristic distribution is integrable and each leaf is symplectic. 

Due to the Leibniz condition there is a correspondence between almost Poisson brackets and bivector 



for /, 5 S C°°{P). We will work indistinguishably with almost Poisson brackets and bivector fields. We 
denote by vrf : T*P TP the map such that /3{J{a)) = n{a,/3). More generally, if (j) £ n''{P), then 
Tr'^{4>) is a fe-vector field defined on 1-forms ai,...,Qfc by 7r''((/))(ai, a^) = (— l)'^(/)(7rS(ai), 7rS(afc)). 

Note that the characteristic distribution of vr is the image of vr" and the hamiltonian vector field 
of / G C°°{P) is Xf = —7r'^{df). The 3-vector field [7r,7r], where [•, •] is the Schouten bracket (see e.g. 
|13|). measures the failure of the Jacobi identity of {•, •} through the relation 



{f9,h} = f{g,h} + {f,h}g, 



for f,g,hGC^{P). 



Xf{g) = {g,f}, 



fields vr G r(/\^ TP) given by 



T^{df,dg) = {f,g} 



- [vr, vr] {df, dg, dh) = {/, {5, h]} + {5, {h, /}} + {h, {/, g}} 



(2.6) 




[a,/3]^ = ^,rti{a)/3 - £-n^(p)Ot - d{7r{a, (3)) = ■£^t!(a)/3 - K«ii3)da. 



(2.7) 
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The bracket [•, ■]tj is M-bilinear, skew-symmetric and satisfies the Leibniz identity. In general, vr" is not 
necessarily bracket preserving; instead (see e.g. O Section 2.2]), 

7r»([a,/3]J = [7rtt(a),7r«(/3)]-^W;3[vr,^], for a, /? G [^^(P). (2.8) 

If the bivector field vr is Poisson, i.e., ['/r,7r] = 0, then there is Lie algebroid structure induced on T*P 
with bracket [•, -Jtt and anchor map vr^ : T*P ^ TP, see e.g. |13] . 

We will use the following formula that characterizes the Schouten bracket of a bivector field. 

Lemma 2.1. If (P, vr) is an almost Poisson manifold, then for f,g,h G C°°{P), 

^[7r,7T]{df,dg,dh) = cyd. [dh{7:\[df ,dg]^)) + dh{[K\df),'n\dg)])\ . (2.9) 

Proof. It is a direct consequence of ()2.6p and the fact that the bracket defined in (|2.7p satisfies \df, dg]-,^ = 
d{f,g}. □ 

A regular distribution F C TP and a nondegenerate 2-section Vtp ^ T{/\^ F*) define uniquely a 
bivector field vr on P by the relation that, for a G T*P and X ^ F, 

\x^F = a\F TT\a) = -X, (2.10) 

where a\F denotes the point-wise restriction of a to F. Observe that, in this case, F is the characteristic 
distribution of the bivector vr. Moreover, any bivector field vr with regular characteristic distribution 
arises in this way. In particular, if F is a regular distribution on P and 17 G r2^(P) is a 2-form such 
that Vt\F is nondegenerate then the pair [F, 0) defines a bivector field on P. Note that in terms of the 
bracket {■, •} associated to vr, (|2.1U|) is written, for /, 5 G C°°(P), as 

{/,5} = f)i.(vrtt(d/),^tt(d5)). 

Given a nondegenerate 2-form (7 on P, (|2.1U|) defines an associated bivector field vr for F = TP. In 
fact, there is a one-to-one correspondence between nondegenerate 2-forms and nondegenerate bivectors 
(i.e., those for which the map tt'^ : T*P — )• TP is an isomorphism). The Jacobiator in this case satisfies 

-dn{ni{df),ni{dg), vr« {dh)) = {/, {g, h}} + {g, {h, /}} + {h, {/, g}}, 

for f,g,h G C°°{P), which is equivalent to 

^[vr,vr]=vr«(df^). (2.11) 
Gauge transformations of bivector fields 

It will be convenient to consider a special equivalence relation between bivector fields defined by the 
so-called gauge transformations [28] , 
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Let vr be a bivector field on P and i? be a 2-form on P such that the endomorphism (Id + o vr") 
of T*P is invertiblci, where B^ -.TP ^T*P IS given by B\X) = ixB. The gauge transformation of 
vr by the 2-form B is the bivector field iTg on P such that 

7ri = ttK (Id + B^ o Jy\ (2.12) 

Two bivectors vr and vr are called gauge related if there exists a 2-form B defining a gauge transfor- 
mation from TT to vf. 

Remark 2.2. 

(i) Gauge related bivectors have the same characteristic distribution. In particular, if vr is a bivector 
field with an integrable characteristic distribution, then any gauge related bivector has the same 
associated foliation, even though the leafwise almost symplectic structures may change. Also 
note that gauge related brackets have the same Casimirs. 

(ii) If vr is a regular bivector field defined by a distribution F and a nondegenerate 2-section i}p then 
a gauge transformation of vr by the 2-form B is the bivector field Tig defined by the distribution 
F and the 2-section Qp — B\f. Note that the invertibility of the endomorphism (Id -|- i?'' o vr^) is 
equivalent to — B\f being nondegenerate. 

(Hi) If TT and vf are bivector fields defined by nondegenerate 2-forms ft and Q (via (|2.10p with F = TP) 
then they are automatically gauge related by the 2-form B = Q — Q. 

(iv) Let TT and vr^ be gauge related bivector fields by a 2-form B. If tt is Poisson and B is closed then 
ttb is Poisson, [28] . 



Example 2.3. On with coordinates {x,y, z,pi,p2), consider the following bivector field: 

d d d d , d d , „n 

ox opi ay op2 opi op2 

where a, b are non-zero functions, and the 2-form B = bdx A dy. The endomorphism Id + B^ o vr" is 
invertible since, written in the basis {dx,dy,dz,dpi,dp2), B^ o vr" is an upper triangular matrix. In 
order to compute the gauge transformation of tt by B, denoted by iig, we use equation (j2.12p to see 
that 

TT^idx) = 4 o {dx + B^o J{dx)) = Triidx) 

Doing the same for the rest of the variables we obtain that '/r|((ipi) = TT^{abdy + dpi) and 7r|((ip2) = 
7r'^{— bdx + dp2), while on the other elements of the basis vr and iTg coincide. Therefore we obtain that 
d d d d 

TTg = a— A 7- h 7^ A TT — . Similarly, the gauge transformation of vr by —B is 

ox opi Oy op2 

d d d d ^ d d , ,s 

^-fl = «ir ^ + TT ^ 7^ 2a6— A— . 2.14 

ox opi oy op2 opi op2 

The three bivectors vr, iTg and 7r_s are gauge related and z is a Casimir for all of them. Finally, observe 
that = ^ A — cib-^ A is not gauge related with vr since x is a Casimir only for vf . 



gauge transformation of a bivector field can be defined for any 2-form but, if tliis invertibility condition is not 
satisfied, the result is an almost Dirac structure 1281. 
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2.2 Nonholonomic systems and the nonholonomic bracket 



Let us consider a nonholonomic system determined by a lagrangian L : TQ ^ M of the form 

C = K + U, 

where k is a Riemannian metric (called the kinetic energy metric) on Q and U € C°°{Q) is the potential 
energy, and a nonintegrable distribution D on Q (i.e., a non-involutive subbundle D C TQ), which 
describes the permitted velocities of the system. We denote the annihilator of D by D°. 

Note that the Legendre transform Leg : TQ — ?> T*Q is a global diffeomorphism since Leg = k'^, 
where k'^ : TQ — > T*Q is given by {k'^{v),w) = k,{v, w) for v,w G TQ. 

Let n € C°°(r*Q) be the hamiltonian function given hy V. = C o Leg"^ : T*Q R. If, locally, 
D° = spanje" : a = 1, k.} where e'* are 1-forms on Q, then the equations of motion can be written, 
in coordinates {q^,Pi), as a first order system on the cotangent bundle T*Q given by 

.,■ dH . d'H. , n • - -, -\ 

Wi' ^^ = "9^ + ^'^'^' ^ = l,---,n, (2.15) 

where Xa,a = 1,. . . ,k, are functions which are uniquely determined by the fact that the constraints 
are satisfied. Then the constraint equations become 

ef{q)jf^=0, a = l,...,k. (2.16) 

OPi 

The geometry underlying a nonholonomic system allows one to write the equations of motion in an 
intrinsic way. Let A4 := Leg(-D) C T*Q be the constraint submanifold. Since the Legendre transform 
is linear on the fibers, is a vector subbundle of T*Q. We denote by r : — >■ Q the restriction to 
M. of the canonical projection tq : T*Q ^ Q. 

If Qq is the canonical 2-form on T*Q, we denote by the 2-form on M defined by '■= i'*^Q 
where l : Ai —?■ T*Q is the natural inclusion. 

The constraints are intrinsically written as a regular, non-integrable distribution C on A4 given by 

C^ = {vC TmM : Tt{v) G Z),(^)}. (2.17) 

A fundamental result given in [H Section 5] is that the point-wise restriction of 0,^ to C, denoted by 
O^lc, is nondegenerate. 

The nonholonomic vector field G X(A^) is the vector field uniquely defined by 

ix,A\c = dTij^lc, (2.18) 

where Hm is the restriction to M of the hamiltonian function H, i.e., = ''*^ : — ?■ M . The 
integral curves of X^h are solutions of the system (|2.15|) . 

We define the nonholonomic bivector field ir^i, to be the bivector associated to the distribution C 
and the 2-form as in (j2.10p (see [29 \ 124 ^ [T9]). Then (A^jVTnh) is an almost Poisson manifold. The 
distribution C is the characteristic distribution of the bivector tt^i, and since C is always not integrable, 
TTnh is never a Poisson structure. The nonholonomic bivector describes the dynamics in the sense that 
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So, in hamiltonian form, a nonholonomic system is described by the triple (A^, vr^h, T^jvi), defined 
from the lagrangian L and the distribution D. 

We call (A4,7rniJ the almost Poisson manifold associated to the nonholonomic system. 

In [2] we saw that it is worth to explore different bivectors tt on A4 such that = —'^'^{d'HM)- In 
this case we say that vf also describes the dynamics. 

In order to generate new bivector fields on A4 we will use gauge transformations of the classical 
nonholonomic bivector n^i-, by 2-forms (as in Section 12. ip . If we want to generate bivector fields 
describing the dynamics, we need the following definition given in [21 Sec. 4. 2], 

Definition 2.4. Let (P, tt) be an almost Poisson manifold with a distinguished Hamiltonian function 
H € C°°{P). Given a 2-form B on P such that (Id + B'" o n'^) is invertible, the gauge transformation 
of TT associated to the 2-form B is said to be a dynamical gauge transformation if 

1x^3 = 0, 

where Xh is the hamiltonian vector field associated to H. 

Thus, any bivector iTg that is dynamically gauge related to vr satisfies 7rg{dH) = ■K'^{dH). We are 
interested in using Definition 12.41 to perform a dynamical gauge transformation of the almost Poisson 
structure vTnh, where the distinguished function is Tij^. 

Remark 2.5. In [21 Proposition 11], it was shown that if we consider a semi-basic 2-form B with 
respect to the fiber bundle r : — > Q ( i.e., ixB = if Tt{X) = 0) then (Id + B^ o vrL) is invertible. 

o 

3 The Jacobiator of the nonholonomic bracket 

Let (A^,7rnh) be the almost Poisson manifold associated to a nonholonomic system. In this section we 
present a coordinate- free formula for the Jacobiator of the classical nonholonomic bivector tt^^, or any 
bivector gauge equivalent to it. This section will not use the hamiltonian function T-Lm- 

Splittings of TM. adapted to the constraints 

Recall that the characteristic distribution of tt^y, is the nonintegrable distribution C defined in ()2.17p . 
Since C is a regular distribution on A4, it is possible to choose a complement W of C on such that, 
for each m G M, 

TmM=Cm®Wm- (3.19) 

Consider the projections P^ : TM. — C and Pw : TM. — > W associated to the decomposition (|3.19p . 
Since Pyv '■ TM — )■ W can be seen as a W- valued 1-form, following [6l [25] , we define the W- valued 
2-form Kyv on M given by 

Ky,{X, Y) = -P^{[Pc{X), Pc{Y)]) for X,Y e X{M). (3.20) 

It is straightforward to check that Kw is (7°°(A4)-bilinear. Note also that Kw = if and only if C 
is involutive (see also [25l Section 2.5]). 
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Proposition 3.1. The W-valued 2-form Kw is semi-basic with respect to the bundle projection r : 
M^Q,i.e., ixKw = i/ Tt(X) = 0. 

Proof. Let us consider X,Y & X(7W) such that X is r-related with and Y E r(C) is r-projectable. 
Since Ker Tr C C then [Pc{X), Pc{Y)] = [X,Y] is r-related with 0. Therefore, [Pc{X), Pc{Y)] G T{C) 
and thus, from the expression ()3.20p we obtain that 'K-^^{X,Y) = 0. This implies that i^Kw = since 
r-projectable vector fields generate TAi at each point. □ 



Gauge transformations adapted to the splitting 

From ()2.12p we observe that a gauge transformation of 7r„ij by a 2-form B is completely determined by 
the point-wise restriction of -B to C (i.e., by the section B\c of /\^C*). That is, if two 2-forms B and 
B are such that {B — B)\c = 0, then vrg and TTg agree. 

Note that for any 2-form B, one can find another 2-form B which agrees with B when restricted 
to C, so that TTg = TTg, and that satisfies the additional condition 

izB = for ah Z G r(W). (3.21) 

So from the point of view of gauge transformations, there is no loss in generality in assuming that 
()3.21|) holds. This condition will be assumed in the sequel to simplify some formulas. 



The Jacobiator of 7r„h and gauge related bivectors 



Now we state the theorem which characterizes the Jacobiator of any bivector -Kg gauge related to n^h- 
Even though the most interesting statement for this Theorem occurs when a Lie group is acting, we 
will first state it in the more general setting without the consideration of a symmetry. 

Let us consider the almost Poisson manifold (7W,7rnh) associated to a nonholonomic system, as in 
Section [221 



Theorem 3.2. Let ttb be a bivector on M gauge related to n^i-, by the 2-form B. Let W be a complement 
of C and assume that B verifies (|3.21|) . Then for a,l3,j G T*A4 we have 

^[7rs,7rs](a,/3,7) = cycl. Jl^(Kw(vr|(a), 7rfs(/3)), 7ri(7)) - 7(Kw(vrfj(a), 7r|(/3)) ; 

+ dS(4(a),4(/3),4(7)). 

Proof. Since this formula is C°°(A^)-linear on a, /3, 7, it is sufficient to check it on exact forms. Applying 
formula ([2^91) we obtain, for f,g,h e C°°{M), 



1 



'TTB,7rB]{df,dg,dh) = cycl. 7Tl{[df,dg]nB){h) + [7Tl{df),7ri{dg)]{h) 



cycl. 



-7rUdh){{f,g}B) + dh{Pci[nUdf),7rUdgm 
+ dhiP^i[7riidf),7ri{dg)])) 



(3.22) 
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where we are using the fact that Id = + Pw- If we call ilj^ '■= ~ B, then by Remark I2.2f zi) and 
()2.10p we have i j /,,h\^M\c = —dh\c and thus 

dh{Pc{[T:l{df)A{dg)])) = -^M{T^l{dh),Pc{[7:l{df)A{dg)])) 

= [4(4f), 4(^^5)1) + f^^(4(d^),^'w([4(d/), 4(^9)])). 

Thus, ()3.22p is equivalent to 

]:['KB,-nB\{df,dg,dh) = cycl. [-7r|(d/i) (o^(7r|(d/), 7r|(d5())) + $7^([7r|(4f), 7r|((i5)], 7r|((i/i)) 

^ L V / ^3 23) 

-o^(Pvv([4(#), 4(^^5)1), 4 (d/i)) -Kvv(4(d/),4(d5))(M , 



where the last term is the definition of Kw given in (j3.20p . 

Note that the cyclic sum of the first two terms of ()3.23p is just — (iO^(4(4f)) 4(c^5)) 4('^^)) 
which is equal to di3(4(4f)) 4(t^9)) 4(c^^))- By (|3.2ip we also have that i 



1 



^M- Finally, again using (j3.20p . we obtain from (|3.23p 



Pv^{[^i{df),J^{dg)\)^^ 



-[TrB,TTB]idf,dg,dh) = dB{7rl{df),TTl{dg),Trl{dh)) 



+ cycl 



(Kw (4 (df) , 4 (dg)), 4 (dh)) - Kw (4 (df) , 4 {dg)m 



□ 



As a particular case of the previous Theorem, we obtain a formula for the Jacobiator of the non- 
holonomic bivector tt^y^. 

Corollary 3.3. The Jacobiator of the classical nonholononiic bivector TTn^ is written, for a choice of 
W, as 

J[vrnh,vr„h](a,/3,7) = cycl rJl^(Kw(4h(a), 4h(/3)), 4h(7)) - 7 ('Kw(4h(a), 4h(/3)) 



In \25\ Section 2.5] , Marsden and Koon presented a local version of the above formula for an explicit 
choice of W that is induced by the local coordinates. Also, in agreement with [25] . we see that ir^i, is 
Poisson if and only if Kw = 0. 



4 The Jacobiator of reduced nonholonomic brackets 

In this section we use Theorem 13. 2l to obtain a formula for the Jacobiator of any bivector field ttb gauge 
related to TTn^ when the nonholonomic system admits symmetries. This formula will be in terms of 
a curvature-like object and the momentum associated to the cotangent lift of the group action. This 
will lead to a formula for the Jacobiator of reduced brackets, and it will provide information about the 
integr ability of their characteristic distributions. 
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4.1 Nonholonomic systems with symmetries 

Let G be a Lie group acting on Q freely and properly. Suppose that G is a symmetry of a given 
nonholonomic system, that is, the lifted action on TQ leaves the lagrangian C and the constraints D 
invariant. Then, M is an invariant submanifold of T*Q by the cotangent lift of the action. Therefore 
we restrict the G- action on T*Q to M.: 

(f : G X M ^ M, 

so CxC"^) = ^T'ginT-)- As a consequence of the G-invariance of the constraints and lagrangian, the 
nonholonomic bivector vr^h and the hamiltonian function are invariant by the G-action ip. 

Throughout this section, we will not use explicitly the hamiltonian function instead we will 
focus on the G-invariant distribution D on Q and the G- invariant almost Poisson manifold (MjTr^i,). 

The dimension assumption 

At each m € Al, let us denote by Vm ■= ^m(Orb(?Ti)) the tangent to the orbit at m of the G-action on 
A4. We will always assume that the following condition, known as the dimension assumption, holds: 

Trr^M =Cm + Vm, for all m G 7W. (4.24) 

Let us denote by S the distribution on given, at each m € by 

Sm = Vmr\C^. (4.25) 

It is clear from ()4.24p that S has constant rank. If, at each q ^ Q, we denote by := Tq{Oih{q)) the 
tangent to the orbit associated to the G-action on Q, (|4.24|) is equivalent to the dimension assumption 
considered in [6]: 

TgQ = Dg + Vg, QGQ, 

by the definition of C and since, at each m G A^, Vm and V^(m) are isomorphic. Analogously, we can 
define the distribution 5 on Q given hy Sq := Dg n Vg for each q € Q. 

If 5 = {0}, the system is called Chaplygin and the dimension assumption is just 

TM = C®V. (4.26) 

From now on, we will always assume that the symmetry Lie group acts freely and properly on Q 
so that the dimension assumption is satisfied. 

The momentum map 

Recall that the submanifold l : M ^ T*Q is invariant under the cotangent lift of a G-action by 
symmetries. Let us denote = t'*^Q and Qm = i'*®Q, where fig and &q are the canonical 2- form 
and the Liouville 1-form on T*Q, respectively. We let : — ?> g* be the restriction to A4 of the 
canonical momentum map of the lifted action on T*Q: 

(Jim),!]) = i^^(„,)G^(?n), for r/ G g. (4.27) 
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Note that = d{J',ri). More generally, we will need to consider vector fields rjjii associated with 

elements t] € T{M. x g) = C°°{Ai,g) which are not constant (i.e., ??A4(m) := {rim)M{^))j this case, 
we have 

ir,j^nj^ = {dj,r]). (4.28) 
Splittings of TA4 adapted to the symmetries and constraints 

Let (7V4,7rnh) be the G-invariant almost Poisson manifold associated to a nonholonomic system with 
symmetries satisfying the dimension assumption (|4.24|) . Let W be a G-invariant distribution on A4 
such that, at each m € A1, 

T,r,M = Cn,®W,n and WmQVm. (4.29) 

Note that Vm = Sm © Wm- We refer to such W as a G -invariant vertical complement of C in TA4. 

Since the projection Pyy, : TA4 W associated to the decomposition (j4.29p is an equivariant map 
(i.e., Pw(Tipg{vm)) = Tipg{Pw{vm))), it induces a g-valued 1-form Aw by 

Aw{Vm) = C Pw{Vm) = S.Mi'm), for Vm G TmM. (4.30) 

Observe that Ayv ■ TM. g is, in general, not surjective; note also that it satisfies Ayv{TLpg{vra)) = 
AdgAw{vm) for g e G. 

Remark 4.1. A choice of a G-invariant vertical complement of D in TQ, 

TQ = D®W and W QV, (4.31) 

induces a G-invariant W satisfying (|4.29p . Indeed, since Trlv : V ^ y is an isomorphism, we can take 
W := (rrlv)-^(T^). 

Note that such G-invariant vertical complements W always exist: for example, one can take W 
to be y n S^, where is the orthogonal complement of S with respect to (the G-invariant) kinetic 
energy metric k, as proposed in [6]. In this case, Ay^, = r*^"^'", where is the g- valued 1-form on 
Q used in [6] to define the nonholonomic connection. However, in this paper, we will allow different 
choices of VF; in fact, we will consider many examples in which W is not defined in this way. o 

The W-curvature 

A G-invariant W as in (|4.29p and the resulting map A^j ■ TA4 — )> g (|4.3Up induce a g-valued 2-form on 
Ai, that might be interpreted as the associated "curvature". More precisely: 

Definition 4.2. The W-curvature is the g-valued 2-form on Ai given by 

ICyv{X,Y) := dAwiPciX),Pc{Y)) for X,Y e TM, (4.32) 
where Pc ■ TM — ?> C is the projection associated to the decomposition (j4.29p . 

Remark 4.3. In the case of a Chaplygin system, where C is a complement of the vertical space V, the 
only choice of W is W = V, hence A\> = ^ is the principal connection on M with horizontal space 
given by C. Therefore /Cy = is the classical curvature. o 
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Proposition 4.4. For any G-invariant vertical complement W of C in TA4, the W -curvature /Cw 
verifies the following: 

(i) }C^{X,Y) = -AU[X,Y]) forX,YG T{C). 

{ii) /Cw is ad-equivariant, i.e., £^^/Cw = ad^lCw, for C S 5- 

Proof (i): U X,Y e r(C) then 

JC^{X,Y) = dA^{X,Y) = X{A^{Y)) - Y{AUX)) - A^{[X,Y]) = -A^{[X,Y]). 

(ii): For g G G, we have 

(^*/Cw(X,y) = dAUPcTipg{X),PcTipg{Y)) = -AU[T^gPc{X),TipgPc{Y)]) 

= -Ay,{T^g[Pc{X),Pc{Y)]) = -AdgAU[Pc{X),Pc{Y)]) = AdgJC^{X,Y), 

so the claim in (i) follows. □ 

As a consequence of item (i) and ()4.30p . we see that for X,Y G X(7W) we have that 

Kw(X,y) = (/Cw(^,n)^, (4.33) 

for Kyv in (|3.20p . Hence, /Cvv = if and only if C is involutive. 

Finally, we consider a natural 3-form associated with the momentum map and the W-curvature. 
Let X,Y,Z be vector fields on M. Then /Cw(X,y) G T{M x g) and dJ{Z) € r{M x q*). We define 
dj A /Cw to be the 3-form on Ai given by 

dj A /Cw(^, Y, Z) = cyclic [{dJ{Z),JC^{X, Y))] (4.34) 

where (•, •) is the pairing between q* and g. 

4.2 The Jacobiator for the nonholonomic brackets in the presence of symmetries 

As a consequence of Theorem 13.21 we obtain the following Jacobiator formula: 

Theorem 4.5. Let (A^,7rnh) be the G-invariant almost Poisson manifold associated to a nonholonomic 
system with symmetries satisfying the dimension assumption (|4.24p . Suppose that -Kg is a bivector 
field on M. gauge related to vTnh by a 2- form B satisfying ()3.2ip . and let W be a G-invariant vertical 
complement ofC. Then 

]^[kb,t^b] = -7ri{dB + dJ AlC^)-ij^^, (4.35) 

where : Ai ^ Q* is the canonical momentum map, /Cw is the W-curvature, and ^t^^ is the trivector 
V'^^(a,/3,7) =cycl [7 (^w(4(a), 4(/3))) ] , for a, e T*M. 
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Proof. Let X,Y,Z G r(C). Note that ICyv{X,Y) E T{A4 x q) is not necessarily a constant section, so 
using ([OS]) and ([03|) . we get 

J^^(Kw(X, y), Z) = J7^((/Cw(^, Y))^ , Z) = {dJ{Z) , JCUX, Y)). 

From Theorem 13.21 and using ()4.34p . the Jacobiator of vr^ can be written as 

^[7rB,TTB]{a,/3,j) = (dS + A /Cw) (vri(a),7r|(/3),7r|(7)) - cyd. 7 (^/Cw(vr|(a), 7r|(/3))^ 
which completes the proof. □ 

Note that the lack of integrability of C can be seen directly from Theorem 14.51 Since C is regular, 
it is integrable if and only if it is involutive. By (|2.8|) and (|4.35|) . one sees that C fails to be involutive 
due to the presence of the trivector ipT^g (which vanishes if and only if /Cw = 0). 

Corollary 4.6. For a nonholonomic system with a symmetry group satisfying the dimension assump- 
tion ()4.24p . and for a G-invariant vertical complement W ofC, the Jacobiator of the nonholonomic 
bivector n^h is given by 

1 a 



4.3 The Jacobiator of reduced brackets 

For a nonholonomic system with symmetries satisfying the dimension assumption, we will now study 
the Jacobiator of reduced brackets on M./G. 



Reduced brackets 



Consider the G-invariant almost Poisson manifold vTnh). Note that if a bivector vrg is gauge related 
to TTnh by a G-invariant 2-form then it is also G-invariant (see [21 Proposition 13]). In this case, the 
orbit projection p : M ^ M/G induces almost Poisson bivector fields ttJ'^'Ji and tTj^^ on Ai/G, where 

TpoT:l{p*a) = «:,)«(«) and TpoJ^{p*a) = (vrf ,)»(«), for a G ^\M/G). (4.36) 

Equivalently, let {•, - js be the almost Poisson bracket on G°°(A^) associated to the invariant bivector 
ttb- The induced almost Poisson bracket {•, -j^^ on functions on M./G is given by 

{fop,gop}B{m) = {f,g}ZM^))^ for /,5 G G^iM/G),m € M. (4.37) 

Let W be a G-invariant vertical complement of C. From Theorem 14.51 we have 

Corollary 4.7. Let ttb be a bivector field on A4 gauge related to 7r„h by a G-invariant 2-form B 
satisfying condition (j3.2ip . Then 

<A = -4{dB + djA /Cw) o p*. 
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Proof. If TTj^j is p- related to Tig then, from ()4.36p . their Schouten brackets are p-related as weh. That 
is, 

[^fcd. Trfd] = Tp[ttb,ttb] o P*- 

□ 

As observed in |2], vr^^JJi and vr^^ can have fundamentally different properties even though tt^i-, and 
ttb are gauge related. For instance, in general 7r°J^ and nf^^ have different characteristic distributions 
(and thus are not gauge related). In the special case that the 2- form B is basic (with respect to the 
orbit projection p : ^A ^ A4/G), then the reduced brackets n^^^ and vr^^^ are gauge related: 

Proposition 4.8. Let P be a manifold equipped with an action of a Lie group G that is free and proper. 
Let IT and ttb be G-invariant bivector fields on P, and denote by vr^ed and vr^^ the induced bivectors on 
P/G. If TT and ttb are gauge related by a basic 2-form B (with respect to p : P ^ PjG), then vr^cd and 
'^fcd ^''^ gauge related by the 2-form B^^^ where B^^^ is the 2-form on P/G such that p*B^^^ = B. 

Proof. Since (Id + B'" o vr") is invertible on T*P then, using the relations ()4.36p and the fact that B is 
basic, we obtain that (Id+S^^^jOvrL) is invertible on T*{P/G). By (l2T^ we have that vrH = 7r|o(Id+5''o 
7r«) which is equivalent to ir^op* = 4 op*(Id+-B^td°7rL). Therefore, = °{^^ + Btd°T^ii)- □ 

The goal of the remainder of this section is to find conditions guaranteeing that reduced almost 
Poisson brackets on Ai /G have integrable characteristic distributions. A class of almost Poisson brack- 
ets that has this integrability property (not necessarily satisfying the Jacobi identity) are the so-called 
twisted Poisson brackets. The appearance of these brackets in nonholonomic systems was initially 
observed in 

Reduced brackets and tw^isted Poisson structures 

A bivector field vr on a manifold P is called a (j)-twisted Poisson j28j if (/> is a closed 3- form on P such 
that 

i[^,7r]=7r«(</>). (4.38) 
If {•, •} is the bracket defined by the (/(-twisted Poisson structure vr, then (j4.38p becomes 
{/, {g, h}] + {g, {h, /}} + {h, {/, g}} - <p{Xf,Xg,Xh) = 0, 

for f,g,heG^{P) and = {•,/}. 

Let TT be a 0- twisted Poisson bivector on the manifold P and [•, -JyT be the bracket on T*P given by 
()2.7p . Note that vr^ does not preserve this bracket. However, using (|2.8p and (|4.38p we obtain 

[7r»(a),7r«(/3)] = vr" ([a, + i^«(„)^^«(^)</>) , 

for 1-forms a, (3 on P. This induces the following modification of the bracket ()2.7p : 
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which makes {T*P, [•, •]</,, vr") into a Lie algebroid [28] (see also [3 Sec. 2.2]). As a consequence, if vr is 
twisted Poisson, its characteristic distribution 7r''(T*P) is an integrable distribution (generally singular) 
since it is the image of the anchor map of a Lie algebroid. Moreover, each leaf lq : O ^ P of the 
corresponding foliation of P is endowed with a non-degenerate 2-form that is not necessarily closed: 
if vr is i;^>- twisted, then L*Q(j) = d^Q. Note also that if vr is twisted and cj)' is a closed 3-form whose 
pullback to each leaf coincides with that of </), then vr is also (/)'-twisted. 

Example 4.9. 

(i) If vr is the bivector field on P defined by a nondegenerate 2-form 0, then vr is ^-twisted and 
(/> = d^, see (|2lT]) . 

(ii) Any bivector field vr on a manifold P with an integrable regular characteristic distribution is 
twisted Poisson, even though there is no canonical choice of closed 3-form (/>, [2]. 

[in) Let vr be a (/)-twisted Poisson bivector field on P. If vr^ is a bivector field obtained from vr by a 
gauge transformation by a 2-form then vr^ is (0 — (ii?)-twisted. In particular, if vr is Poisson, 
then vr^ is (— (ii?)-twisted. 

Coming back to our setting, observe that the almost Poisson manifold {M^tt^^) associated with a 
nonholonomic system is never twisted Poisson since its characteristic distribution C is not integrable. 
Moreover, since all bivector fields gauge related to vr^h have C as characteristic distribution, none of 
them is twisted Poisson. 

Suppose that G is a symmetry group for the nonholonomic system satisfying the dimension as- 
sumption ()4.24p . Consider the corresponding G-invariant almost Poisson manifold (A^jVr^iJ. The next 
result gives sufficient conditions for reduced brackets to be twisted Poisson based on Corollary 14.71 

Corollary 4.10. Let W be a G-invariant complement of C in TA4 such that the 3-form dj A /Cyy is 
closed. Let B he a G-invariant 2-form on M satisfying condition (|3.21|) and defining a bivector field 
TTg gauge related to vr^^. // dB + dj A /Cw is basic with respect to p : M ^ M/G then the induced 
reduced bivector field vr^^ is (p-twisted Poisson for the 3-form cj) on Ai/G defined by the condition 

p*cl) = -{dB + dJ MCy,). 

In particular, if we now consider the G-invariant hamiltonian T-Lm : A1 — > M and if the 2-form B in 
Corollar V 14.101 defines a dynamical gauge transformation (as in Def. l2.4p . then the reduced nonholonomic 
vector field X^^^^ can be written as 

where T-L^^^ : Al/G — M is the reduced hamiltonian function ( i.e, T-Lm = P*'H^cd)- So, in this case, 
the reduced dynamics is described by a twisted Poisson bracket. Note that we are not assuming any 
regularity condition on vr^^^. 

4.4 A special class of symmetries 

We now show how many of the results discussed so far simplify once we impose an additional condition 
on the G-invariant complement W of C 
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We say that W satisfies the vertical-symmetry condition if there is a Lie subalgebra 0w ^ Q such 
that, for aU m € Al, 

W„ = {eMM|eG0w}. (4.39) 

Observe that we only need to require gyv to be a vector subspace since in this case it automatically 
satisfies that [g, gw] C flvv by the G-invariance of W. 

The vertical-symmetry condition is a restrictive condition on the symmetries, as it implies the 
existence of a subgroup Gw of G for which the system is Chaplygin. But, as we will see, many 
non-Chaplygin systems satisfy this condition, including the vertical rolling disk (Section 17. ip . the 
nonholonomic particle (Sections 14.151 I5.3|) . the snakeboard (Section I7.2p and the rigid bodies with 
generalized constraints (Section 17. Sp ). 

Remark 4.11. Note that systems satisfying the vertical-symmetry condition can be seen as Chaplygin 
systems with additional symmetries, so one can in principle treat them in two steps: first taking care 
of the Chaplygin symmetries, then treating the extra ones (along the lines of e.g. |181 l27j). Comparing 
this viewpoint with the direct, one-step approach is the subject of a separate work in preparation [T]. 
o 

The next result shows some consequences of (|4.39p . 

Proposition 4.12. For any G-invariant vertical complement W of C in TM satisfing the vertical- 
symmetry condition ()4.39p . the W -curvature /Cw verifies 

(i) JC^{X, Y) = dA^{X, Y), for all X G T{C),Y G X{M). 

{ii) dlCyv{X,Y,Z) = 0, for all X,Y,Z er{C). 

(Hi) dj A = d{J , JCyv) ■ 

Proof {i): U X,Y £ r(C), then ]Cyv{X,Y) = dAyv{X,Y) by the definition of /Cw If X G T{C) and 
Y G r(W), then /Cw(^,>") = 0. In order to check that dAwiX,Y) = 0, we write Y = fii])^ where 
fi G C°°(A1) and if G 0w Then, using that [X, r/^] G T{C) (because C is G-invariant) and that 
•^w{Vm) = ry*, we get 

dA^{X,Y){m) = dA^{X,f,r]i,) = {X{rj') - A^{[X,7f^])) = 0. 

(ii) : Let X,Y,Z G r(C). Then using item (i) we obtain 

d/Cw(X, Y, Z) = cycl[Z{lC^{X, Y)) - }C^{[X, Y],Z)] = cycl[-Z{Ay,{[X, Y])) - dAy,{[X, F], Z) ] 
= cycl[-Z{Ay,{[X, Y])) + Z{Ay,{[X, Y])) + AU[[X, Y], Z\) ] = 0. 

(iii) : It follows from ()4.34p and d/Cw|c =0. □ 

We can simplify Theorem 14.51 and Corollarv 14. 71 when W is chosen such that condition ()4.39p holds. 

Theorem 4.13. Let (A1, TTnh) ^'e the G-invariant almost Poisson manifold associated to a nonholonomic 
system with symmetries satisfying the dimension assumption ()4.24p . Let W be a G-invariant vertical 
complement of C satisfying the vertical- symmetry condition (|4.39p . Then: 
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(i) the Jacobiator of a bivector field iTg on A4 gauge related to vr^jh by a 2-form B satisfying condition 
()3.2ip is given by 

1 J, 

-[7rs,7rs] = -7r%{dB + d{J,K.w)) - V'tts- 
(ii) //, in addition, B is G-invariant, then the Jacobiator of the reduced bivector field nf^^ is 

Ik,, <J = -^UdB + d{J, /Cw)) o p*. 

Corollary 4.14. Under the same hypotheses of Theorem \^.l'J\ (ii), if B is such that dB + d^J^K-w) 
is basic, then nf^^ is a (p-twisted Poisson bivector on M/G where cj) is the 3-form on M/G such that 

P*(l)=-{dB + d{J,ICyy)). 



Example 4.15 (The nonholonomic particle). Consider a particle in Q = subjected only to the 
constraint e = 0, where 

e = dz — ydx, 

with lagrangian C = ^{x^ + + i^). The system admits symmetries given by the action of the Lie 
group on Q such that the vertical space isV = span { ^ j ^ } • The non-integrable distribution D on 
Q is given by -D = span |^ + 2/^) ^| and we choose the complement W oi D to he W = span {^}- 
Observe that W is not equal to S-^ H V. In canonical coordinates, the constraint manifold is 



M = {{x,y,z;px,Py,Pz) ■ Pz = VPx}- 



From here we see that 



C = span<^ — + y— — ^ and W = span <^ — ^ . 4.40 
[ ox oz oy opx opy J [oz ) 

Note that the vertical-symmetry condition ()4.39p is satisfied, for 

= span{(0, 1)} C = R\ (4.41) 

Remark 4.16. Note that the subbundle Qs := {{m,(^) : & <Sm} C A4 x q m this case is given 

by 0s I m = span{(l,y)}, for m £ M; so the system does not admit horizontal symmetries in the sense 
of [6]. So, the vertical-symmetry condition may hold without horizontal symmetries. o 

The projection to W adapted to the decomposition (|4.40p is given by Pw{vm) = ^i'^m)-§^, for 
Vm G T„iM, where we consider {dx,dy,dpx,dpy,e} a local basis of T*^A. Therefore Aw ■ TM. — > 
is given by Ay^ = (0,e) (see ()4.30p ). Since 0vv is a vector space, the W-curvature is /Cw = (0,(ie) = 
{^,dx A dy). To compute the momentum map J : M. ^ q* we use (|4.27p for the canonical 1-form 
©Ai = Pxdx + pydy + ypxdz and (1,?/)ai = ^ + and (0,1)^1 = Therefore, in canonical 
coordinates, ^7 = ((1 + y'^)Px,yPx)- The 2-form (j7, /Cyv) is given by 

{J,]Cw) = yPxdxAdy. (4.42) 
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In this case, we can explicitly check that dj NK-y^ = di^J ^ A^w)- By Theorem 14 . 5 1 or 13^. 1 3 1 we have the 
ingredients to write the Jacobiator for Tr^h (and for gauge related bivector fields tts). More interesting 
is analyzing the Jacobiator of the reduced bivector field vr^^Jj on the reduced manifold M. /G ~ (with 
local coordinates iy,Px,Py))- Observe that (J7,/Cw) is not basic (not even d{J',ICw) is basic) but 

d{J, ]C^){7TUp*dy),7TUp*dp^),7TUp*dpy)) = 0. (4.43) 

By Theorem l4.13l fii). we conclude that [tt^^^^, tt^^\] = 0, i.e., the reduced nonholonomic bracket is Poisson. 
Equation (j4.43|) can be easily checked from the fact that 



dpy ' " ' y"^ + 1 \dx dz J + 1 dpy ' " dy y"^ + 1 dpx 

More examples where d{J',ICyv) is basic, or where we need a dynamical gauge in order to obtain 
d{{J', K-Yj) + B) basic, are worked out in Section [71 



4.5 Chaplygin systems 

As we now see. Theorem 14.131 and Corollary 14. 141 have direct applications to Chaplygin systems. 

For a nonholonomic system with Chaplygin symmetries (see Section 14. ip , the splitting (|4.29|) be- 
comes 

TM = C®V, {AAA) 

since S = {0}. That is, the only possible choice for W is exactly the vertical space V, which clearly 
satisfies the vertical-symmetry condition ()4.39p . In this case /Cyv = \s the usual curvature of the 
connection associated to the splitting (|4.44p . 

Recall that for a Chaplygin system, the reduced bivector ■k'^^^ is defined by a nondegenerate 2-form 
ri"^*"^ on M./G. A direct consequence of Corollarv 14. 141 is the following result of [1]: 

Corollary 4.17. The 2-form 0."^]^^ satisfies 

dnf:, = -d{J,JC)..„ (4.45) 
where (J',/C)rod "is the 2-form on A4/G such that p* {{J , K) red) = {J-,^)- 



Proof. Starting with Corollary 14.141 we set B = and note that {J', K) is a basic 2-form, since it is 
G-invariant and semi-basic with respect to the orbit projection p : 7W — )• M./G. It then follows that 

which is equivalent to (|4.45p by (|2.1ip . 

□ 



Remark 4.18. The Theorem stated by Bates-Sniatycki in [1] was written in terms of the horizontal 
lifts, i.e., dO.^'^^iX.Y^Z) = -d{J,IC){X^ , Z^) where X,Y,Z e X{M/G) and , Z^ G r(C) 

are their horizontal lifts to C o 
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One also has a more general statement including gauge transformations. Since W = V, condition 
(j3.2ip reduces to requiring that B is semi-basic with respect to A4 M /G: 

izB = forallZGV. (4.46) 

So if B is G-invariant, it must be basic, and there is a 2-form B,.^^ on M/G such that p*B^^^ = B. 

If TTg is gauge related to TTn^ by a basic 2-form B, then (by Proposition 14. 8|) their reductions tt^^^ 
and vr^^j are gauge related by -Bred- Since tt"^]^ is defined by a nondegerenate 2-form 0,'^^^, it follows 
that vr^^j is defined by a nondegenerate 2-form ^f^^, which is related to O"^*"^ by = ri"^''^ — B^^^ (see 
Remark 12. 2l f in)). 

Taking into account gauge transformations, Corollary 14. 141 savs that 

l[^L,<J = -i^LfidB,^, + d{J, /C),e,) 
which, by (j2.1ip . is simply the assertion that 

dnf^, = -d{{j,ic)^,, + B,^,). 

5 Poisson structures associated with some nonholonomic systems 
with symmetries 

For a given a nonholonomic system with an action by symmetries satisfying the dimension assump- 
tion, this section shows how the choice of a G-invariant vertical complement W of C satisfying the 
vertical-symmetry condition endows the reduced manifold M/G with a genuine Poisson structure. We 
will also dicuss conditions under which the reduced nonholonomic dynamics is described by a gauge 
transformation of this Poisson bracket. 

5.1 Poisson structures on reduced spaces 

Let (A^,7rnh) be the almost Poisson manifold associated with a nonholonomic system equipped with a 
symmetry group G such that the dimension assumption is satisfied. We start by observing 

Proposition 5.1. For any G-invariant vertical complement W ofC, the 2-form {J', ICyv) on M. satisfies 
the following: 

(i) it is semi-basic with respect to the bundle projection t : M ^ Q, 
{ii) it is G-invariant. 

Proof. Statement (?) relies on (|4.33p and on the fact that Kw is semi-basic, which was proved in 
Proposition 13.11 

To prove {ii), recall that /Cvv is ad-equivariant (see Prop. By the equivariance property of the 

canonical momentum map J : AA ^ Q* obtain 

^e.M^'^ '^w) = (-^5a4^ '^w) + {J , £(m^w) = {-ad*^J,}Cw) + {J ,ad^lCw) = 0. 

□ 
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As a consequence of the previous Proposition and Remark l2.5l if B := —{J', /Cw), the endomorphism 
°'^th) of T*A4 is invertible. Therefore, the gauge transformation of ir^^^ by the 2-form — {J^ , /Cw) 
is well defined and produces a new bivector field on ^A that we denote by ttj^ • 

Equivalently, we can describe ttjic by noticing (see Remark 12. 2( ii)) that the point-wise restriction 
of the 2-form 

to C is nondegenerate. Then ttjx is the bivector determined by the distribution C and the section ^jx\c 
(as in (|2.10p ). Note also that ttj/c may not describe the dynamics since A^vv) usually does not 

vanish. 

Proposition 5.2. Assume that the G-invariant vertical complement W of C satisfies the vertical- 
symmetry condition ()4.39p . Then ttj^: is a G-invariant bivector field satisfying 



and its reduction tt^ defines a Poisson structure on M/G. 

Proof. The G-invariance of ttj^: follows from the G-invariance of {J',fCyv), see Prop. IS.lf u). The 
formula for [ttjic ^t^jk ] follows from Theorem I4.13r i). once we set B = —{J',ICyv)- The fact that vr;^ 
is a Poisson structure is an immediate consequence of Corollary 14.141 □ 

To simplify our notation, we will denote the Poisson structure on A4/G defined in Prop. [5^2] bv A, 

A := TT^, , (5.47) 

assuming that the choice of W is clear from the context. 

Example 5.3 (The nonholonomic particle). Recall Example 14.151 The nonholonomic bivector tt^i^ is 
given by 

y"^ -\- 1 \dx ^ dz J dpx dy dpy 1 + dpx dpy 

Observe the similarity with the bivector in (j2.13p . The gauge transformation of vr^h by the 2-form 
-(J,/Cw) given in (I02]l . is 

1 / d d \ ^ d d ^ d 2 ypx d ^ d 
^'^'^ y'^-\-l\dx ^ dz J dpx dy dpy l + y'^dpx dpy 

The orbit projection p : Ad ^ Ai/G is given by p{x , y , z; px , Py) = iy,Px,Py) and thus the Poisson 
bivector field on M./G is 

dy dpy l + y^ opx opy 



Example 5.4 (Chaplygin systems). In the case of Chaplygin symmetries, as recalled in Section | 
the reduced bivector tt^^^ on M/G is defined by a 2-form 1^"^''^, and the 2-form {J ^IC) is basic with 
respect to the orbit projection p : M. ^ M./G. Hence there is a 2-form {J,1C)^^^ on M./G so that 
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/C)rcd = (i7,A^)- By Prop. 1121 the reduced bivector Aon ^4/G is obtained from vr^^'jj by a gauge 
transformation by —{J'jIC),.^^, which is equivalent to saying that A is defined by the nondegenerate 
2-form 

n':^, + {J,JC)^^,. (5.49) 

The assertion that A is a Poisson structure is the same as saying that the 2-form ()5.49p is closed (i.e., 
symplectic), which is in agreement with (|4.45p . 

5.2 A metric-independent Poisson manifold 

In this section we only use part of the data of a nonholonomic system with symmetries, since the 
lagrangian will play no role. 

Consider a (free and proper) G-action on the configuration manifold Q preserving the constraint 
distribution D, and let be a G-invariant vertical complement of D in TQ, that is, TQ = D (BW 
and W CV. 

Analogously to the condition introduced in Section 14.41 we will say that W satisfies the vertical- 
symmetry condition if there is a Lie subalgebra Qw of Q such that, for all q Q, 

Wg = {vQ{q)\v<^Qw}. (5.50) 

Since W is G-invariant, the same holds for its annihilator W° C T*Q with respect to the lifted 
action to T*Q. We will now show that W°/G inherits a Poisson structure Aq when W satisfies the 
vertical-symmetry condition. Later, in Section 15.31 we will relate the Poisson manifold {W°/G, Aq) 
with the reduced Poisson structures of Prop. 15.21 

Let Vo C T{W°) be the vertical space with respect to the G-action on W° (the restriction of the 
cotangent lift), and denote by tq : W° Q the canonical projection (the restriction of tq : T*Q Q). 
As in Section 12.21 let us consider the distribution Cq on W° defined by 

Co := {v G T{W°) : Ttq{v) G D}. (5.51) 

Observe that Cq is G-invariant, since D is G-invariant. 

Consider Ovi^o := tg^g, where lq : W° ^ T*Q is the natural inclusion. 

Lemma 5.5. The restriction ^lw°\co ^-^ nondegenerate. 

Proof. The lemma is a direct consequence of the result in [H Sec. 5] recalled in Section 12.21 if we 
consider any metric kq for which D and W are orthogonal, then we have A4 = Kq{D) = W°, and the 
induced constraint distribution on is Cq. □ 

Therefore, as seen in Section \2.2\ the distribution Cq, given in ()5.5ip . and the 2-form Qy^^ induce a 
G-invariant bivector field vro on W° defined by the relation 

ix^w-\co = a\co ^ 7r^(") = -^ at each q G r*(T^°). (5.52) 
We denote the reduction of ttq to W°/G by Aq. 
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Then, from Remark 14. 11 we know that the decomposition TQ = D@W induces a sphtting T{W°) = 
Co ©Wo, where Wo := (rro|vo)~"^(^)- this case, since W is G-invariant then Wo is G-invariant and 
also note that if W verifies the vertical-symmetry condition, then Wo verifies it as well. 

Proposition 5.6. If W is a G-invariant vertical complement of D satisfying the vertical- symmetry 
condition, then Kq, the reduction of ttq on W° to W°/G, is a Poisson structure. 

Proof. By Theorem 14. 13| we know that 

1 



2 



ko,vro] = -■Kl{d{Jo,ICwo)) -i^- 



•TO ) 



where j7o : W° — ?> g* is restriction of the canonical momentum map on T*Q to W° and /Cwq is 
the Wo-curvature with respect to the decomposition T{W°) = Cq (B Wo- We are going to show that 
(i7o)^wo) = Oi for which it suffices to check that (,7o; ^wo) Ico = 0- Consider (local) basis {Xi,Za} 
of TQ adapted to the decomposition D © W. If {q) are local coordinates on Q, then we denote by 
{q;Pi,Pa) the coordinates on T*Q associated to the dual basis {X^,Z"'} ofT*Q. The submanifold W° 
is represented by {q;pi,0). In order to compute the 2-form (j7o,/Cwo) observe that for X,Y G r(Co), 
Pvvo([^,^]) = Z''{[X,Y]).Za = -dZ''{X,Y).Za. If is a local basis of sections of g^,. x Q such that 
{^^)q = Za then the coordinates of K,y^,,{X,Y) are (K.y^^,{X,Y))a = dZ''{X,Y). On the other hand, 
note that the Liouville 1-form on T*Q can be written as 

0Ql(9;p„p.)=P*^*+Pa^". (5.53) 

(Here we use the same notation for 1-forms on Q and their pullbacks to T*Q.) So (JojC^) = i-oPa- 
Hence (i7o,^wo)lco = i^oPa)dZ"'\cQ = 0, since iQPa = 0. Therefore, applying Theorem \A.13( ii). it is 
clear that the reduction of ttq is Poisson. □ 

We now use the results in the Appendix to give an alternative characterization of the Poisson 
structure Ao on W°/G. Consider the 2-form Qyyo on W° and the distribution Wo as in ()5.5ip . 

Lemma 5.7. IfW satisfies the vertical- symmetry condition, then Keril.w° = Wo- 

Proof. For v € Wo|m> ^ ^ W°, we have that Qw°{'v){m) = m{TTo{v)) = 0, since Ttq{v) € W. Hence 
Wo C Ker0,yo. To see that Wo C KevQ^^, it is sufficient to check that i^^^Jlwo = for all ?7 € g^, 
where r^vvo is the infinitesimal generator corresponding to rj. Then 

using the G-invariance of Qw° ■ As remarked in Lemma 15.51 \co is nondegenerate, hence Wo = 
Ker . □ 

It follows from Lemma lA.lf i) in the Appendix that the presymplectic submanifold {W°,^}w°) of 
{T*Q,Qq) reduces to a Poisson structure on W°/G. The fact that this Poisson structure coincides 
with Ao of Prop. 15.61 is a direct consequence of Lemma lA.41 

Proposition 5.8. Let W be a G-invariant vertical complement of D in TQ satisfing the vertical- 
symmetry condition. Then the Poisson structure on W°/G obtained by reduction of (as in 
Lemma \A.H i) ) coincides with the Poisson structure Aq of Prop. \5.6[ 
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Example 5.9 (Chaplygin systems). In the case of Chaplygin symmetries, W = V so the Poisson 
structure Aq is defined on V°/G. By Prop. 15. 8|, Lemma 15.71 and Remark IA.2I in the Appendix, Aq 
corresponds to a symplectic form on V° /G, uniquely determined by the property that 

p*ilo = = i^o^Q- (5.54) 

If J : T*Q — > 0* is the canonical momentum map for the lifted action on T*Q, then V° = J~^(0) 
and it is well-known that the corresponding Marsden-Weinstein reduction is the cotangent bundle 
(r*(Q/G),OQ/G), see e.g. [26]. With the identification V°/G = J-^{0)/G = T*{Q/G), it follows from 
(|5.54p that coincides with the canonical symplectic form 0,q/g, so the Poisson structure Aq is just 
the canonical one on T*{Q/G). 

5.3 The equivalence 

In this section, we compare the Poisson structures A and Aq, obtained in Propositions 15.21 and 15.61 We 
recall the set-up defined by a nonholonomic system with symmetries. 

Let D be the G-invariant constraint distribution on the configuration manifold Q, and let be a 
G-invariant vertical complement of D. This defines an invariant bivector field ttq on the submanifold 
W° C T*Q as in (|5.52|) . On the other hand, by using the (G-invariant) kinetic energy metric k, we 
consider the submanifold A4 C T*Q equipped with the G-invariant vertical complement 

W := {Tt\v)'\W) (5.55) 

of C (induced by W as in Remark 14. ip . On Ai , we have the invariant bivector field ttjjc , introduced in 
Section [5. 1[ Note that if W satisfies the vertical-symmetry condition, then W also does (and gw = Qw)- 

The next result relates the almost Poisson manifolds {A4,ttjx ) (which depends on the metric k) 
and (VF°,7ro) (which is metric-independent). 

Proposition 5.10. For a nonholonomic system with symmetries satisfying the dimension assumption, 
let W be a G-invariant vertical complement of D in TQ, and consider the almost Poisson manifolds 
{J^,TTj^ ) and {W°,-Ko). Then there exists a G-equivariant diffeomorphism ^' : — )• W° preserving 
the almost Poisson structures. 

As a consequence, ifW satisfies the vertical- symmetry condition, so that {A4/G,A) and {W°/G,Aq) 
are Poisson manifolds, ^ induces a Poisson diffeomorphism between them: 

iM,7Tj^) 1 > (VF°, vro) (5.56) 

{M/G, A) {W°/G, Ao). 

Proof. Let us consider a local basis {Xi, Za} of TQ adapted to the decomposition TQ = D (B W, such 
that {Xi} is orthonormal with respect to the kinetic energy metric k. Let {X^, Z"} be the dual (local) 
basis of T*Q. This basis induce local coordinates {q;Pi,Pa) on T*Q, where (q) are local coordinates on 
Q. The Liouville 1-form Qq on T*Q is written as Qq = piX^ + paZ"" (see ()5.53p ). so 

J^Q = A dpi - p^dX' + Z'^ A dpa - PadZ". (5.57) 
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The local coordinate description of A4 = k?{D) is {{q^Pi^ f^iaPi)}^ where nia = K{Xi,Za), since 
/{"(Xj) = X* + K{Xi, Za)Z"'. So the pull back of the canonical 2- form Q,q to A4 is 

nM=X' A dpi - Pi ■ dX' + A d{KiaPi) - (KiaPi) ■ dZ"". 

Observe that, locally, C = spanjXj, ^ + Kjq^-} and W = spanjZa}. 

Let us compute the 2-form {J', /Cw) (as done in the proof of Prop. 15. 6p in these coordinates. We 
take {£,""} a local basis of sections of Qw x Q such that {(,'^)q = Za- For X,Y r(C), the coordinates 
of JCyv{X,Y) are (/Cw(X,y))a = dZ''{X,Y). Therefore we obtain 

{J,IC^)\c = i^iaPi)-dZ''\c, (5.58) 

and, since Z"|c = 0, 

^jk\c = i^M + (J,/Cw))|c = {X' A dpi-pi ■ dX')\c. (5.59) 

Consider the metric kq on Q given at each X,Y G TQ by 

Ko{X,Y) = k{Pd{X),Pd{Y)) + k{Pw{X),Pw{Y)), (5.60) 

where Pd '■ TQ D and Pw '■ TQ W are the projections associated to the decomposition TQ = 
D(BW. The spaces D and W are orthogonal with respect to this metric. Moreover, the basis {Xj, Za} is 
KQ-orthonormal. Then the coordinate description of W° = Kq{D) is {{q',Pi,0)}, and the decomposition 
T{W°) = Co ® Wo is given, locally, by Cq = spanjXj, ^} and Wo = spanjZa}. On the other hand, 

= t*nQ = X' A dpi - Pi ■ dX\ (5.61) 

We have to see that there exists a diffeomorphism ^' : — > W° preserving the almost Poisson 
structures ttj^ and ttq. For that, it suffices to check that 

T*(C)=Co and Vtj^\c = ^*{VLwo)\c. (5.62) 

The map 

^> := {kq\d) o : M ^ W° (5.63) 

is a diffeomorphism that, by the definition of k, satisfies ^'(X* + KiaZ"") = (Ko|D)(Xt) = X*. That 
is, the global diffeomorphism ^ : M ^ W°, written in local coordinates, maps {q;Pi, HiaPi) on Jli 
into {q',Pi,0) on W°. A direct computation shows that the conditions in ()5.62p hold, so ^ preserves 
almost Poisson structures. It follows from (|5.63|) that ^ is G-equivariant so that the induced map 
M-jG ^ W° /G also preserves the reduced brackets. 

□ 

Example 5.11. Let us keep the notation of Example 14.151 and 15. 3i Consider the local coordinates 
{px,Py,Pz) associated to the local basis {dx,dy,e} of T*Q. Then, {dx,dy} generates W° so the asso- 
ciated local coordinates in W° are {px,Py,0)- Therefore, 

r^w" = dx A dpx + dy A dpy = dx A dpx + dy A dpy 
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( d d d d d ] 

and Co = span < — — h 2/7:-, t:-, — , — r where (px,Pv,Pz) are the canonical coordinates on T*Q. 
[ ox oz ay opx opy J 

The isomorphism ^ : ^A ^ W° of Proposition 15.101 written in canonical coordinates, is 
^{Px,Py,yPz) = (^Y^^'-PyO^ • We can check that ^'*Ow°|co = {'i^ + y'^)dx Adp^+dy Adpy + lyp^dx A 
dy = Ic- 
on the other hand, the Poisson bivector field Aq is computed using the presymplectic 2-form 0,^^° 
and ()A.78p and thus we obtain that Aq = ^ A The Poisson bivector A in (|5.48p is isomorphic to 
Aq via the isomorphism '^.^^ : M/G ^ W°/G. 



We now combine Prop. 15.101 with Examples 15.41 and 15.91 to recover the following result from |22j 
about Chaplygin systems (see section 

Corollary 5.12. On a Chaplygin system, we have an identification M/G = T*{Q/G) in such a way 
that 

where I^q/g is the canonical symplectic 2-form on T*{Q/G). 

Proof. By Example 15.41 the Poisson structure A on Ai/G corresponds to the symplectic form ^l'^^^ + 
/C),.<,d, while Example 15.91 shows that W°/G = V°/G, equipped with the Poisson structure Aq, is 
identified with the canonical Poisson structure on T*{Q/G). By Prop. 15.101 once we identify M/G 
with V°/G = T*{Q/G), A coincides with Aq, which is equivalent to Q'^^^ + {J, JC),,^ = Qq/g- □ 



6 Symplectic leaves and the reduced dynamics 

Throughout this section, we consider nonholonomic systems with symmetries satisfying the dimension 
assumption. We let W Q TQ be a G-invariant vertical complement of the constraint distribution D 
satisfying the vertical-symmetry condition (|4.39p . and let W C TA4 be as in Remark 14.11 Our aim is 
to relate the Poisson structure ()5.47p with the reduced nonholonomic dynamics. 



6.1 Symplectic leaves 



Let us consider the Poisson structures A on Ad/G (as in ()5.47p ) and Aq on W° /G (as in Propositions [5i6 
and 15. 8p . We start by describing their symplectic leaves. 

Using the description of Aq in Prop. 15.81 as a reduction of Vtyv°i we have the following diagram: 



1-0 



->(r*Q,QQ) 



(6.64) 



PO 



(T^7G,Ao)^--> (T*Q/G,7r.. 



As remarked in Lemma lA. 11 the symplectic leaves of Aq are the (connected components of) the inter- 
sections of the leaves of vr^an with the submanifold W°/G T*Q/G. Moreover, the leaves of Aq sit 
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inside the leaves of tt^^^ as symplectic submanifolds. It is worth noticing that W° /G is not generally 
a Poisson submanifold of T*Q/G (it is a Poisson-Dirac submanifold in the sense of \10\ Sec. 8]). 

More explicitly, denoting hy J : T*Q q* the canonical momentum map for the lifted action on 
T*Q and letting j7o = Jolq : W° g* be its restriction to W°, it directly follows from Cor. lA.3l that the 
symplectic leaves of Aq are given by (connected components of) the quotients J7q~^(//)/G^, for fi G g* . 
The symplectic structures on the leaves are described as follows: recalling that the symplectic leaves of 
TTj-an are given by (connected components of) the Marsden-Weinstein quotients (J~^(/i)/G^, w^), [261 
Sec. 2] the symplectic structure on Jq^{^)/G^ is given by Z^w^, where we keep the same notation 



for the restriction of the the inclusion. 

Using the isomorphism ^ : — ?> W° of Prop. 15.101 and the induced Poisson diffeomorphism 
*rod : (Al/G, A) (VF°/G, Ao), we conclude the following. 

Theorem 6.1. Consider the natural embedding lq o : A4/G ^ T*Q/G. Then: 

(i) The symplectic leaves of the Poisson manifold {A4/G,A) are connected components of the inter- 
sections of leaves o/vr^an with the image of A4/G in T*Q/G, and they sit inside the leaves o/vTean 
as symplectic submanifolds. 

(ii) The symplectic leaves of A are connected components of quotients {Jq o ^)^-'^(/.i)/G^, for n € 
0*, equipped with symplectic forms {Iq o ^^^^)*uj^, where we keep the notation lq o "^^^^ : (^q ° 
^')~"'^(/i)/G^ ^ for the restricted embedding. 

We illustrate this result with the nonholonomic particle, following Examples 14.151 15.31 and 15. Ill 

Example 6.2. Recall that, in canonical coordinates, T*Q/G is given by {y,Px,Py,Pz) and the leaves of 
the Poisson bracket vTcan are the level sets of the functions px and pz. Since W° /G is locally described 
by the coordinates {y,Px,Py), then, as we just saw, the leaves of Aq are the level sets of px. Now, using 

the diffeomorphism ^ : Ai ^ W°, given in canonical coordinates by '^ipx,Pv,yPz) = ( — o iPv-i^] i 

we obtain the symplectic leaves of A: by Theorem 16.11 they are given by the inverse image by ^rcd of 
the leaves of Aq, that is, the level sets of (1 + y'^)px determine the leaves of A. 

6.2 Relation with the nonholonomic momentum map 

In what follows, we analyze the role of the nonholonomic momentum map [B] in the description of the 
symplectic leaves of the Poisson manifold {M/G,A). 

The definition of the nonholonomic momentum map, that we now recall, does not use any choice 
of complement W of D. Following [6], let us consider 5 = C R V C TAi (see ()4.25p ) and define the 
subbundle — t- A1 of the trivial bundle = A4 x g ^ A4 by 



We can similarly consider S = V Ci D C TQ and Qs C QX Q, Qs\g = : S,q{q) £ Sq}. Since the action 
on M. is the lifted action on Q, Qs\m coincides with Qs\g for r(m) = q. 



io:Jo\fi)/G^^J~\fi)/G^ 




(6.65) 
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For an arbitrary section ^ G ^{Qs) C r(g^) = C°°{M-,g), we have 

Hm^m = -kMd'S>M = d{^j^@M) - £^M®M- (6.66) 

The nonholonomic momentum map is the map JT""'" : M ^ Q% given by 

J''\mm = kM{m)®M{m), (6.67) 

for C G flsU- Note that J'* G r(g*) so that, if ^ G T{qs), then = Hm^m G C°°(>t). 

From equations (|6.66p and (|6.67p . we see that, for any ^ G r(g5), we have 

i^^n^=d{J'^\0- £(m®m. (6.68) 

If we now consider the G- invariant hamiltonian function G C°°(A4), then from ()6.68p we observe 
that (^"'',0 € C°°{M) is conserved by the flow of X,,,, if and only if ix^^^^£^j^^@M = 0. 

Remark 6.3. If ^ is a constant section, then £^^0^1 = because of the G-invariance of the Liouville 
1-form. This implies the conservation of the nonholonomic momentum map in the case of horizontal 
symmetries, as it was proven in [6]. o 

The next result relates ^7"'' : — )■ and j7o o ^ : — )• g* (considered in Theorem I6.1f n)) once 
a choice of W is made. Observe that, for each m G Al, the Lie algebra g can be split into 

= 0s|m©0w, (6.69) 

observing that the splitting is simplified by the vertical-symmetry condition, in that gw is independent 
of m. Let 

Pqs ■Q^4-^ i3s (6.70) 
be the projection associated to the splitting (|6.69p . 

Remark 6.4. One can consider a similar splitting to (|6.69p if the vertical-symmetry condition does 
not hold for W, except that the complement of g^ will also vary with m; the projection (|6.70p is also 
defined similarly. o 

Theorem 6.5. Consider the almost Poisson manifold {A4,7rj,^), defined in Section \5.1[ The nonholo- 
nomic momentum map J''^^ : A4 satisfies the following: 

(i) Jo°'^ = o J""^. In other words, 

{Joo^{m),7]) = (J'"'^(m),Pg5(r/)), for all m € M,r] £ q. 

(ii) If W verifies the vertical-symmetry condition ()4.39p then vr^((i(^"'', Pg^ (?]))) = —{Pgs{'n))M, 
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Proof, (i) Let us denote, as usual, lq : W° T*Q and l : M. ^ T*Q the corresponding inclusions 
and let Pd '■ TQ D is the projection associated to decomposition TQ = D ® W. Let v £ TQ and 
m € M., then since ^'(m) G W°, we have that {Lo{"^{m)),v) = (io(^'(?n)), ^^(f )). Therefore, if 77 G g, 
and using that 7]q = TTo{r]w°) we obtain 

(Jo ° 'I'M,??) = {Lo{^{m)),i]Q) = {Lo{^{m)),PD{vQ)) 

Now, by (|5.63p . we have that io{"^{'m)) — i{m) € D° C T*Q since if m = K^(ug) where Vq G Dq then 
(io(^'(m)) - /,(m))|B = (J^ivq) - = by Therefore, 

{io{^{m)),PD{vQ)) = {i{m),PD{vQ)) = {i{m),{P,^v)Q) = ( J-'^M, Pg^ (77)). 

(a) Since Jo '■ W° — )■ g* is the restriction of the canonical momentum map on T*Q to W° 
then at each rj G q, i^^^Q^o = d{Jo,rf). By Lemma [5.71 we have that Kerrivv'" = VVo, and de- 
noting by Pco • T{W°) Co the projection associated to decomposition T{W°) = Co © Wo we 
obtain that \p^^(^^^„)Q,w°\co = d{Jo, 11)100. Using (|5.52|) we obtian that 7rl{d{Jo,r])) = -Pco(??w°)- 
Prop. EIO] asserts that, for each a € T*{W°), T'^{Jj,c{^*{a))) = 7rJ(a), and thus vr?^((i(Jo ° ^,?/)) = 
TT^j^{<i>*d{Jo,r])) = -Pc{7]m) since T^PcM) = ^'coM)- Finally, using that PcM) = (^gs(^))Ai 
item (i), we get that Pg^ (??)» = -(n5(^))Ai- 

□ 

Corollary 6.6. The symplectic leaves of the Poisson manifold {A4/G, A) are (the connected components 
of) iJ^^^r\Llfi)/G^, for fi e g*, where l, : ^ is the natural inclusion. 

In the previous statement, we view /i € g* as a constant section of g^, so that is a (not 
necessarily constant) section of g^. 

Proof. By Theorem 16. 11 the leaves of A are the connected components of (Jo ° ^)^^(^)/G'po for fJ- & 9* ■ 
Using Theorem 16. Sf i) and the fact that l* o P*^ = Idg* , we see that (Jo o = ( J"'')~^(/,*^), for 

/ueg*. □ 

Remark 6.7. Following Remark 16. 4[ we note that part (i) of Theorem 16.51 still holds without the 
vertical-symmetry condition. o 



6.3 Reduced nonholonomic brackets as gauge transformations of Poisson struc- 
tures 

As we already remarked in section 15.11 the Poisson bivector A on /G may not describe the reduced 
dynamics, since ttjic does not necessarily describe the dynamics defined by the hamiltonian Hj^. Still, 
there are situations in which the reduced dynamics is described by a bivector field obtained as a gauge 
transformation of A: 

Theorem 6.8. Suppose that B is a G-invariant 2-form on M satisfying ()3.2ip and defining a dy- 
namical gauge transformation of tt^i-, (as in Def. \2.4\ ). If B + (J,/Cvv) is basic with respect to the orbit 
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projection p : A4 A4/G, then the reduced dynamics is described by a bivector field vTj^^, which is a 
gauge transformation of A by the 2-form B, where B is such that p*B = B + (JT, /Cw). 

Assuming that W verifies the vertical- symmetry condition (j4.39|) . we conclude that 7r"^'jj is a twisted 
Poisson structure that is gauge related to the Poisson bivector A. 

The various bivectors and gauge transformations in Thm. 16.81 are represented by the fohowing 
diagram: 

(M, TT^ ) '''"^^'''^^ > {M, vr^) — ^ (A^, 7r„,) (6.71) 
p 



{M/G, A) 2 , (M/G, n 



^ 1 

red / 



Proof. Proposition l4.8l asserts that the bivector vr^^^ is gauge-related to A by B. If W verifies the vertical- 
symmetry condition ()4.39p . then by Corollarv 14.141 bivector Tr^^^ is {—dB)-twisted Poisson, where B is 
the 2-form on M/G such that p*B = B + {J, /Cw). □ 

The procedure in Thm. [6^ will be illustrated in section [7] (see section [7l3]) . 

Corollary 6.9. Under the assumptions of Theorem 1 6'. for each r] & Q, the function -Pg^ (^)) G 

C°°{^A) is conserved by the flow of X^^,. 



Proof. Let ^ = -Pg5(?7), for rj e Q. Theorem |6^u) asserts that 7r2nc(rf(J'°'', 0) = ~C 

j^. oince TTg IS 

the gauge transformation of ttj^ by i? + {J , JC^), we have vr^ = ir'^g o (Id + {B + {J, /Cw))'' o n^j^ ). 
Therefore, 

-U = 7Tlcid{J''\0) = ^kd{J''\0 - kM^B + (J,/Cw») = 7ri{d{T^\0), 

since B + {J',ICyv) is basic. But on the other hand if B defines a dynamical gauge then 7r|(d'H^) = 
—Xnh- Finally, we obtain that ix„h(c^( 0) = —dUM^^M) = because of the G-invariance I-Lm- D 

We stress that the fact that the reduced dynamics is described by a bivector field tTj^^ gauge 
equivalent to the Poisson structure A has several consequences: for example, tt^^ is a twisted Poisson 
structure, so its charctaeristic distribution is integrable - in fact, its almost symplectic leaves coincide 
with the leaves of A, and the Casimirs of A are conserved quantities of the nonholonomic system. In 
particular, the leaves of irf^^ are described by Thm. 16.11 and Cor. 16.61 The only difference is on the 
leafwise 2-forms: while for A the symplectic form on a leaf [Jq o ^)^'^[p)/G^ = {J''^'^)^^ (i^Ip) / G ^ is 

:= {lq o ^>^^^)*uj^ (see Thm. l6.1f M)). for vr^^ the almost symplectic form on the same leaf is 

where B^ is the pullback of the 2-form B on M/G (in Thm. 16. Sp to the leaf. 

We derive some conclusions about the nonholonomic particle based on the results of this section. 
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Example 6.10 (Nonholonomic particle). Recall from Examples 14.151 ^5.3\ 15. IH 16.21 that the reduced 
nonholonomic bivector field on M /G is given by 



.nh 

red 




d d ypx d d 



vr, 



dy dpy 1 + 1/2 Qp 



and the Poisson bivector field A described in Sec. 15.11 is 




1 + 2/2 dpy ■ 



yPx d d 




As we observed in Example 4.15 the 2-form {J', /Cvv) = UPxdx Ady relating vr^h and ttjjc is not basic and 
thus 7r"pj and A are not necessarly gauge related (see Prop. 14. 8p . In fact, for (a, € g ~ and using 
()4.4ip . we see that the function (JT'"'", Pg^ (a, &)) = (J'"'', a(l, y)) = (1 + y'^)px determines the leaves of 
A (Corollarv 16. 6p since / = (1 + y^)px is a basic function on A4. In other words, / is a Casimir of A. 
However, / is not a Casimir of vr" ^ since it is not conserved by the motion (see [5l Sec.5.6]). Therefore, 
TTj^j is not gauge related with A. 

On the other hand, if there was a dynamical gauge B satisfying the conditions in Theorem 16.81 then 
B + {J',lCyv) would be basic and ix^h^ — 0- Hence, ixnh(i7, ^w) would be semi-basic with respect to 
the orbit projection p : M. ^ M./G. However, using equation ()4.42p it can be checked that ^X-^^^ (iT, ^w) 
is not semi-basic. 

As a final attempt to obtain a reduced bivector field describing the dynamics that is gauge related 
to A we may try to change the vertical complement W . Note, however, that in general the vertical- 
symmetry condition will no longer be satisfied. 

7 Examples 

So far we have illustrated our results using the nonholonomic particle (see Examples 14.151 15.31 15. IH 
16.21 and 16.101 and Chaplygin systems (see Section 14.51 and Examples 15.41 and 15. 9p . This section presents 
other mechanical examples, all of which exhibit symmetries which are not Chaplygin, so their reduced 
dynamics are described by bivector fields (not 2-forms). In all cases, we find suitable choices of comple- 
ments W and dynamical gauge transformations so that the reduced dynamics is described by twisted 
Poisson structures (as e.g. in Corollarv I4.1(JI and Theorem 16. 8p . 

More concretely, in Section 17.11 we verify that the reduced dynamics of the vertical rolling disk is 
directly described by the Poisson structure A. In Section 17.21 we consider the Snakeboard, and show 
that the 2-form (c7, /Cw) in this case is basic and induces a gauge transformation relating the reduced 
bivector field vrj* to the Poisson bivector A. Section 17.31 illustrates the usage of a dynamical gauge 
transformation in Theorem 16.81 for the Chaplygin ball, also leading to a description of the reduced 
dynamics by a bracket given by a gauge transformation of the Poisson bracket A. 

7.1 The vertical rolling disk 

Consider a vertical disk of radius R rolling on a plane without sliding, see e.g. [5l Sec. 1.4 and 5.6]. 
Let (x, y) G represent the contact point between the plane and the disk, Lp and represent the 
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rotation angle of a point in the disk with respect to the vertical axis and the orientation angle of the 
disk, respectively. So, the configuration manifold is given by x S*^ x S*^. The nonsliding constraints 
are given by the 1-forms 

Ex = dx — R cos ipdip and Cy = dy — R sin ipdip. 
The Lagrangian is given by the kinetic energy 

y, V'), {x, y, (f, = y(x^ + y^) + ^(199^ + Jip"^), 

where m is the mass of the disk, I is its moment of inertia about the axis perpendicular to the plane 
containing of the disk, and J is the moment of inertia about an axis in the plane of the disk (both axis 
passing through the disk center). Consider the action of the (direct product) Lie group G = x 
on Q given by 

(a, b;a) X (x, y,(p,tp) i-> {x + a,y + b,(p + a, ^). 

The Lie algebra g is M^, and V = span{^, ^, ^}. 

The distribution D given by the annihilator of €x and €y is given below, together with our choice 
of complement W C V: 

W = span\—, — \ , D = span\ — + Rcosi;— + Rsmi;—, — }. (7.72) 
[ ax ay J [dip ox ay dip J 

Note that W verifies the vertical-symmetry condition (j5.50p since 0,^ = span{(l, 0, 0), (0, 1,0)}. 

The dual basis of T*Q associated to (|7.72|) is given by {e^;, e^, dip, dip} with coordinates {px,Py,Pip,Pip)- 
In these coordinates, A4 is defined as 

M = {{x,y,(p,'ip;px,Py,Py,,p^p) ■■ px = {mR^ + I)~^Rmcos^p^, py = {niR^ + ly''^ Rm cos ^ p^} . 

We define, as usual, W as in ()5.55p . thus obtaining the decomposition TM = C W. 
The canonical momentum map J : M ^ Q* associated to is 

J = ii^*Px, t'*Py,Pip) = {{mR^ + ly^Rmcosipp^ , (mi?^ + Rm cos ip p^ , p^), 

where t : — > T*Q is the inclusion. The map Aw is M^-valued, Ayv = {ex,Ey,0), and thus the W- 
curvature is /Cw = {R sin ip dip Adip , —R cos ip dip A dp , 0) . Finally, we see that {J', /Cw) = 0. Therefore, 
the bivector ttjx. coincides with vTnh and thus the reduced dynamics is described by the Poisson bivector 
A (|5.47p on the manifold M./G (given by coordinates (^,P(^,p^)). 

On the other hand, the submanifold W° of T*Q is locally expressed as W° = {px = Py = 0}, and 
thus the presymplectic 2-form il^o is given by Qw° = dip A dp^ + dp A dp^p. Using Prop. 15.81 on the 
reduced space W°/G we have the Poisson structure Aq = ^ A which is computed using (|A.78p . 
To compute the Poisson bivector A on A4/G we consider the isomorphism ^ : A4 ^ W° given by 

Px 5 Py ) Pip 1 

) = {q;0,0,pp,p^). Thus, by Proposition |5JiH 

A- A A — 

dtp dp^ ' 
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so is clearly a Casimir. 

Next, we analyze this Casimir in the context of conserved quantities. At each m & Ai, = 
span{(i2cos'(/', Rsimp, 1)} (see (|6.65|) ). Then for r/ = (a, 6, c) € g ~ M^, the function 

{J''\PsAv)) = (J"'\c(i?cosV,iisin^,l)) =cp^, 

is conserved during motion by Corollarv l6.9l (in agreement with [5] , where it is shown that if {p^ ,Py,Pip,Pi()) 
are the canonical coordinates on T*Q, then = Cp,p is conserved where C is a non-zero constant). 

7.2 The snakeboard 

The snakeboard is a modified version of the skateboard which is modeled on the configuration manifold 
Q = SE(2) X X S^. The coordinates {x,y,9) € SE(2) represent the position and orientation of the 
center of the board, cj) € is the angle of the momentum wheel relative to the board, and G is the 
angle of the back and front wheels (we are assuming that they coincide). The non-sliding constraints 
are given by the annihilators of the 1-forms 

ex = dx + r cot (j) cos 9d9 and Cy = dy + r cot (p sin 9d6. 

The (direct product) Lie group G = x S"-^ is a symmetry of the system (see [27]) with the action 
on Q given by 

(a, b;a) X {x, y, 9, (p,^;) t-^ {x + a,y + b, 9,(l),il} + a,) 
{ d d (9 1 

and thus V = span < — , — , 7— > . Consider the following basis of D and the vertical complement W 



given by 



id n .9 . n . 9 d d } , „r (8 8] 

D = sp&n i — - r COS 9 cot (p- r smt' cot 0— , — , — ^ and VK = span <^ — , — ^ . 

[ o9 ox oy ocp dip J yox ay ] 

In this case the vertical-symmetry condition (|5.50p is satisfied for g,^ = span{(l, 0, 0), (0, 1, 0)}. Hence, 
it is sufficient to compute the 2- form {J',ICw), for W as in ()5.55p . The lagrangian function is just the 
kinetic energy metric, 

£ = 1 (m{±^ + y^ + r^9^) + 239tP + 2Ji(/.2 + J^-^ 

For q = {x,y,9,(j),'ip) G Q, let us denote {px,Py,Pe,P(i},Pip) the coordinates on T*Q associated to the 
basis {€x,ey,d9,d4>,d'tlj}. Therefore, the manifold M = Leg(Z)) is given in coordinates by 

f mr cos 9 sin^ cp cot cp mr sin 9 sin^ (p cot (p ] 

M = \ px = 5 T • 2 - yPe - P^)' Py = 5 t ■ 2, W -Pi^)) ■ 

y mr^ — J sm (p mr^ — J sm (p ) 

Since Qs = span{(0, 0, 1)}, the momentum map adapted to this basis is, for m E A4, J{m) = 
{t-*Px, i'*Py,Pip), where l : Ai ^ T*Q is the inclusion. The W-curvature is given by /Cw = {dex,dey,0), 
and thus the 2-form {J', /Cw) is 

TTiT^ cot (f) 

{J, /Cw) = 5 , . 2 A ^PS - P^)dG A d(p, 

mr'^ — J sm m 
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which is basic with respect to the orbit projection p : M ^ MjG. By CoroHarv 14.141 the reduced 
bivector field vr^^'Ji is twisted Poisson, and hence has an integrable characteristic distribution. Moreover, 
by Theorem 16.81 7r"^'Ji is gauge related with the Poisson bivector A on Al/G by the 2-form (JT", /Cw)rcd 
(here we are taking i? = 0). In order to obtain the conserved quantities of the system that are related 
with the symmetries, we analyze the Poisson manifold (VF°/G, Aq). The submanifold W° C T*Q is 
described in coordinates by p^; = = 0. Using (|5.57p we see that Vl^ = d9/\dp0 + d(j)f\dp^ + d'il)^.d'p^. 
Its reduction by G gives the Poisson bivector Aq on W°/G (Prop. 15. 6|) described, in local coordinates 
{e,4>,pe,P4„p^), by 

A ^ A ^ A ^ 

" 89 dpe dcj) dp^' 

for which is a Casimir. The Poisson bivector A is isomorphic to Aq via the map ^'red defined in 
(|5.56|) and thus '^*^^p^ = p^ is also a Casimir of A (Theorem 16. ip . Since 7r"jjj and A are gauge related, 
the leaves of vr^^'jj are given by the level sets of p^ and thus the function p^ is a conserved quantity of 
the nonholonomic system, in agreement with [27]. On the other hand, observe that, for € 3 — M'^, 
the function p^ satisfies that p^ = {J'''\ Pgs iv)) = (^"'\ (0,0, 1)), see Corollary ESI 

7.3 Rigid body with nonholonomic constraints 

We now consider the motion of an inhomogeneous sphere whose center of mass coincides with its 
geometric center that rolls without slipping on the plane. 

Following [21 Section 5], we can see this example as a particular case of a rigid body with generalized 
rolling constraints which describes the motion of a rigid body in space that evolves under its own inertia 
and it is subject to constraints that relate the linear velocity of the center of mass with the the angular 
velocity of the body. If x G represents the center of mass and u the angular velocity of the body 
with respect to an inertial frame, then the constraints are written as 

k = rAu, (7.73) 

where r denotes the radius of the sphere and A is a given 3x3 matrix that, following [2], satisfies one 
of the (simplified) conditions: 

/000\ /010\ /010\ 

A = 0, A = i , A = I -1 , or A = i -1 . (7.74) 
\001/ \000/ \001/ 

We can interpret the motion determined by (j7.73p for each matrix A: if yl = 0, (|7.73p represents 
the motion of a free rigid body; If A has rank 1 then the body is allowed to move along the vertical 
axis satisfying the constraints; When the matrix A has rank 2, (|7.73p determines a ball rolling on a 
plane without sliding (the Chaplygin ball), while if the rank of A equals 3 then any two points on the 
configuration space can be joined by a curve satisfying the constraints. 

The configuration manifold is Q = S0(3) x with coordinates (g,x), where g is an orthogonal 
matrix that specifies the orientation of the ball by relating two orthogonal frames, one attached to the 
body and one that is fixed in space, and x represents the position of the center of mass in space. We 
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will assume that the body frame has its origin at the center of mass and is aligned with the principal 
axes of inertia of the body. These frames define the so-called space and body coordinates, respectively. 

Given a motion {g{t); x{t)) G Q, the angular velocity vector in space u and body coordinates fi are 
respectively given by 

u{t) = g{t)9^Ht) n{t) = g~\t)g{t), 

where we use the usual identification ' : so(3) — )• M"^, and satisfy u = gfl. Therefore, the constraints 
can also be written as x = rAgfl. 

Let A and p be the left and right Maurer-Cartan forms, respectively, that can be thought as Re- 
valued 1-forms by means of the identification of the Lie algebra so (3) with by the hat map. Thus 
for a tangent vector Vg E TgS0{3) we have u = p{vg) and 12 = X{vg). Therefore, the constraints can be 
written in terms of a M^-valued 1-form e = (ei, €2, £3), where ei are the 1-forms defining the constraints. 
In local coordinates, e is given by 

e = dx — rAp = dx — rAgX, 

where we are using the relation p = gX. 

The constraint distribution D on TQ is given by 

£) = span|x"^'^' + rA^| = span |x'°'' + r^c/^| , (7.75) 

where X^^ht ^ (Xf^^t^ fright ^ ^ogM^j moving frame of S0(3) dual to p = {pi,p2,P3); similarly, 

X'°'' is the moving frame of S0(3) dual to A. 

The left action of the Lie group G = {(/i, a) € S0(3) x : he^ = 63} on Q is given by 

(/i, a) : {g, x) ^ {hg, /ix + a) G S0(3) x M^. 

This G-action on Q is a symmetry for the nonholonomic system since the lifted action on TQ leaves 
the constraints and the hamiltonian invariant (observe that hA = Ah when ^4 is a matrix of the form 
described in (|7.74p ). In this case, the vertical space V is given by 

l^-spanl-v-X'^^-x^A+^iA AAA! 

[ ox^ ax^ ox^ ox'^ ox'^ J 

where 7 = (71, 72, 73) is the third row of the matrix g € S0(3). Note that S = span |7 • (X'°^* + rAg-^)^. 
The Lie algebra g associated to G is abelian Lie algebra identified with M x M^. 

Let us define the G-invariant vertical complement W oi D on TQ given by 

VF = span|A A A\ 
\ dx^ ' dx'^ ' dx^ / ' 

The distribution W satisfies the vertical-symmetry condition ()5.50p for g,^ = span{?7* := (0;ej),i = 
1,2,3} where are the canonical vectors in M^. 

If we consider the basis {X'°'* + rAg-^, ^} of TQ then {A, e} is its dual basis on T*Q. Let us 
denote by ((7,x;K,p) the coordinates on T*Q adapted to this basis. 



37 



The Poisson structure Aq on W° /G. 

The submanifold W° of T*Q is represented, in coordinates, by (g,x;K,0), and the G-action on W is 
given by (/i, a) : (5, x; K) 1— )■ {hg, hx + a; K). Thus, (7, K) are local coordinates of the reduced manifold 
W°/G -S^x M3. Since 9^° = K • A, then 

Qwo = A, A J^i - K • dA = Ai A Ki + K1X2 A A3 + K2X3 A Ai + K3X1 A A2. 

Since W satisfies the vertical-symmetry condition, the reduction of the presymplectic manifold {W° ,i}w°) 
gives the Poisson bracket Aq on W°/G described in coordinates (7, K) as 

{Ki,Kj}Ag = -eijiKi {'yi,Kj}Ag = -eijm {7i,7j}Ao = 0- (7-76) 



The Poisson structure A on M./G via the isomorphism : M. ^ W° . 

The lagrangian is given by the kinetic energy 

/:((7,x;n,i) = K((n,i),(n,i)) = l(m).n+^||x|p, 

where I is the inertia tensor which is represented as a diagonal 3x3 matrix whose positive entries, 
Ii , I2 ) I3 are the principal moments of inertia. 
As it was computed in [2] and [T7] . 

>1 = {(5,x;K,p) G r*g : \^ = m + mr'^{Ag)'^ {Ag)n, f> = mrAgn}. 

Therefore, since we consider coordinates adapted to the constraints, the isomorphism \I' : — > W° 
is written in coordinates as x; K, p) = ((?,x;K,0). On the reduced manifold M/G, with local 
coordinates (7, K), the Poisson structure A coincides with the formulas given in (|7.76p . for W as in 

Remark 7.1. At this point, we can check whether the reduced nonholonomic vector field belongs to 
the characteristic distribution of A. In fact, since ([2l I17j) 

we see that A^(fi • dK) = X'^l\. Even though f2 • dK is not necessarily a closed 1-form (observe 
that it depends on the matrix A by the relation between K and fi), the fact that X^^^ belongs to 
A»(T*(7W/G)) says that it makes sense to think that the reduced dynamics may be described by a 
gauge transformation of A. o 



The 2-form {J^Ky,). 

The map Ay^ : TM. — ?• g defined in (|4.30p is given, as in the other examples, by Ayj = (0,e). Using 
Definition 14.21 the W -curvature /Cvv is 

/Cwlc = dAvM\c = (0, -rAg ■ dX)\c, 
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where d\ = -(A2 A A3, A3 A Ai, Ai A A2). 

Now, we have to write the momentum map : Ai q* in the appropriate basis. Let us consider 
the local basis {7 • (X'"^* + rAg-^), af^'} of ^ adapted to decomposition V = 5©>V. Consider 

the section G r(g) such that .^g € r(5). Therefore, the decomposition 5 © W induces a basis of 
(local) sections of r(g) given by {^*^, -q^ ^rf , rj^}. On the other hand, since a (z Ad C T*Q can be written 
in local coordinates as Og = K^X + mr{AgQ,)'^€, then {J{aq),£^^) = ag(7 • (X''''' + rAg-^)) = (K,7) 
and {J'{aq),ri^) = ag(^) = mrKgii)"^ A'^]i. If we denote the first and second rows of the matrix g by 
a = (01,02,03) and /3 = {(^i, (^2, then (gfl)'^ = {{a, ft), {f3,ft), (7,^)) and thus the momentum 
map JT" G T{M x q*) written with respect to the basis {C^,??*} is 

J={{K,'r);mr{gnfA^). 
Therefore, the 2-form {J', /Cw) on A4 is given by 

{J, /Cw) Ic = r^m{gftfA^Ag ■ X x X\c (7.77) 

Now, it is possible to compute the Jacobiator of tt^i-, and tt^^^ by Theorem 14.131 In order to see if 
the reduced bivector ir'^^^ has an integrable characteristic distribution (or if there is a dynamical gauge 
transformation such that vr^^ has an integrable characteristic distribution) it is necessary to analyze 
whether the 2-from {J, /Cvv) is basic with respect to the orbit projection p : M. ^ M./G and as a result 
we will be able to explain the dynamical gauge transformations chosen in |17^ [2]. 

If A = d. In this case, (J",/Cvv) = and thus [vTnhjVTnh] = since iI'-k^^ = (observe that /Cw = 0). This 
is coherent with the fact that A = Q describes the free rigid body. 

If A has rank 1. In this case, A = A3 with A3 = 63 eg', where 63 is the third canonical vector in 'M? . 
Then, using also that (^7 = 7 x A \n\ [2] the 2-form in ()7.77p becomes 

{J, l^w)\c = mr^igOfA-sg • A x A|c = r^m{^, 0)7 • A x X\c = r'^m{-r, 0)7 • dj x dj\c, 

which is basic. From CoroUarv 14.141 the reduction of TTni, induces a bivector vr^^'jj on A4/G that is 
(— /Cw)red)-twisted where (c7, /Cw)red = r^m(7,n)7 ■ d"/ x d-y € ^}'^{M/G) (which explains the 
result in [21 Theorem 6]). In particular, it has an integrable characteristic distribution. Moreover, 
following Theorem 16.81 t^is bivector Tr"^\ is a gauge transformation of A (given in (|7.76p ) by the 2-form 

// A has rank 2. Since A = Id — A3 where Id is the identity matrix, then 

/Cw) Ic = (r2m(n, A X A) - r'^mla. n)j ■ dj x d-f) \c. 

In this case, d{J,lCw) is not basic. Note that one can conclude more: there is no closed 3-form (p 
on Ai/G making the reduced bivector field vrj'^'jj twisted Poisson. In fact, if such 4> existed then, by 
Corollary we would have ^(^7, /Cw)|c = —p*4'\c^ which would imply that \xd{J ,lCyv) = for any 
X G r('5). However one can verify that this is not the case. As a consequence, since vr^g^ is regular 
([2]), it follows that its characteristic distribution is not integrable, in agreement with |17j . 

Let us consider the 2-form B on Ai given by i? = — r^m(ri, A x A), which is the non-basic part of 
(i7, /Cw)- Observe that B satisfies the conditions of Definition [23] and thus B defines a dynamical gauge 
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for the bivector 7r„h and the hamiltonian Tij^ = • [21 eq.(48)]. Moreover, since B is G-invariant, 
it defines an invariant bivector field tts that describes the dynamics. FoUowing Corollary 14.141 the 
reduction of tts gives a (— dSrcd)-twisted Poisson bivector field vrf where p*B^c^ = {J^ /Cw) + -B, i.e., 

Bred = —r^m{'y, fl)f ■ d'y x ^7. 

So the reduced bivector vr^^^ has an integrable characteristic distribution and, moreover, vr^^^ is gauge 
related to the Poisson bivector field A, given in ()7.76p . by the 2- from B. Again, here we are recovering 
the results in [21 Proposition 2 and Theorem 5]. 

If A has rank 3. Finally A^A = Id, then (^7, /Cw)|c = r'^m{fl,X x A)|c. In this case, /Cw) is not 
basic but it defines a dynamical gauge, that is ixnh(^, ATw) = 0. Thus, the gauge transformation of 
TTnh by i? = — (i7, /Cw) gives a bivector -Kg describing the dynamics. Moreover, -Kg coincides with ttjic. 
So, the Poisson bracket A on M/G describes the reduced dynamics (this result recovers [21 Proposition 

!])• 

Conserved quantities. 

For r/ = (/i; (a,6,c)) G g, we have that (c7°'', -Pg^ (r/)) = {J'^'',^^^) = /i(K,7). Since (K,7) is a basic 
function, then it is a Casimir of A according to Corollary 16.61 Since, for any rank of A, the reduced 
dynamics is described by a bivector field nf^^ that is gauge related with A, (K, 7) is also a Casimir 
for vr^^j. Thus, as Corollarv 16.91 asserts. (K,7) is a conserved quantity of the nonholonomic system (in 
agreement with |17j). 

The nonholonomic bracket. 

Observe that it is not necessary to compute the nonholonomic bracket Tr^h in order to obtain its 
Jacobiator and draw conclusions about the integrability of the characteristic distribution associated to 
the reduced nonholonomic bracket. Note that the 2-section ^^mIc and the distribution C defining vrn^ 
are computed to be C = span jX'"^* + rAg-^, ^} and 

n^\c = '^*n,,.\c - {J,JCy,)\c = {Xi/\Ki-K- dX) \c - (J,/Cw)|c. 
From here we can recover the formula for the nonholonomic bracket given in [21 Prop. 15]. 

A Appendix 

Let {R, 0) be a symplectic manifold equipped with a free and proper action of a Lie group G preserving 
f]. There is an induced Poisson structure on R/G by reduction, that we denote by ttu/g- We use the 
notation Oji/c for its symplectic leaves and pp.:R^ R/G for the quotient map. 

Let L : P ^ Q he a G-invariant submanifold of R and ftp := L*il.. We denote by Vp Q TP the 
distribution tangent to the orbits of the restricted action on P, and hy pp : P P/G the quotient 
map. 

Lemma A.l. Suppose that Kerfip has constant rank and Kerfip C Vp. Then: 
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(i) There is an induced Poisson structure ttp/q on P/G. 

{a) Each symplectic leafOp/o of-Kp/c is a symplectic submanifold of a leafOp/c of-Kp/c, given by a 
connected component of the intersection Op/c H {P/G). 

Note that we have an induced embedding 1 : P/G ^ R/G. Condition [ii) says that P/G is a 
Poisson-Dirac submanifold of R/G, in the sense of [lOl Sec. 8] (see also [8j Sec. 6]). We will refer to 
the Poisson structure t^p/g in {i) as the reduction of the presymplectic form Vtp. 

Proof, (i): Consider C°°{P)adm. ■= {/ G C°°{P)\df (Kern p) = 0}. The condition that KerQp has 
constant rank implies that G°°{P)adm is a Poisson algebra with Poisson bracket given by 

{f,g} = -dg{X), 

where X is any vector field on P such that ix^p = df (such X always exists). If C°°{P)^ denotes the 
space of G-invariant functions on P, then condition Kerilp C Vp implies that G°^{P)^ C C°°{P)adm- 
We now check that if /,5 G C°°(P)'^, then {f,g} G G°°(P)'^. For any 7 € G, using that 7*J]p = VLp 
and 7*/ = /, we see that if \x^p = df, then 7*X — X € Keiilp. Hence dg{X) = dg{'~f^:X) when 
g G C°°{P)'^. From that we directly verify that 'y*{f,g} = -'y*{dg{X)) = -dg{X) = {f,g}. So the 
Poisson bracket on C°°{P)adm restricts to C°°{P)^ , and there is an induced Poisson structure vTp/c on 
P/G since p*p : G'^{P/G) G°°{P)^ is an isomorphism. More explicitly, 

TTp/aid) = Tpp{Z) if and only if Z G X{P) verifies that iz^p = —Pp{oi). (A. 78) 

(ii): By [lOl Sec. 8], showing that the inclusion T : {P/G,-Kp/g) ^ {R/G, TTp/c) is a Poisson-Dirac 
submanifold (which is equivalent to the statement in (ii)) is the same as showing that the subbundle 
graph(7r|/c) = {(7r|/c(a), a) | a G T*{P/G)} of T{P/G) © T*{P/G) equals the subbundle 

{{Y, fa) G T{P/G) © T*{P/G) \ a G T*{R/G) s.t. 4/G(a) = ^^(^)}- (A.79) 

Since both subbundles have the same rank (equal to the dimension of P/G, as both are Dirac structures, 
see [8]), it suffices to verify that the subbundle ()A.79p is contained in graph(7rp/(j). 

Let TT be the bivector field corresponding to the 2-form Q on R. Given {Y,L*a) in the subbundle 
()A.79p at Pp{x) = Pp{x), for x E P (Z R, we know that T/9fj(7r''(p*a)) = vr{j./Q(a) = TT{Y), where Ppa is 
taken at x. For Z G T^P such that TppZ = Y, we see that 7r'^{ppa) — Tl{Z) is vertical (i.e., tangent to 
a G-orbit), and hence must be tangent to P since P is G-invariant. It follows that ^^{ppOL) = Ti{X), 
for X G T^P, and Tl{Y) = Tpp{Tt{X)) = Tl{Tpp{X)), i.e., Y = Tpp{X). Finally note that 

ix^p = L*i^t(p*a)^ = -'*P*pOL = -p*p{fa), 
which says that Y = Tpp{X) = 'Kp/Q{i*a), that is, {Y,L*a) G graph(7rp/g). 

□ 

Remark A. 2. In the special case when KerOp = Vp, then the reduced Poisson ttp/g is defined by 
a nondegenerate 2-form ilp/a- Indeed, we can check that ttp/g is nondegenerate by noticing that, if 
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TTp/Q{a) = 0, then there is a vector field Z G Vp = Kerfip such that iz^p = —pp{a) = 0, hence a = 0; 
that is, vTp/Q is an isomorphism. Note also that ilp/c is uniquely determined by the condition that 

p*pVtp/G = ^p. (A.80) 



To verify this last claim, recall that 7rpyg(a) = Tpp{Z), where Z is any vector field on P such that 
iz^p = —pp{a), and that 1tpp{z)^p/g = — It follows that Pp{ij'ppi^z)^p/G) = —Ppa = iz^p, which 
shows that ()A.80P holds. o 

Suppose now that the G-action on R is hamiltonian, with momentum map J7r : i? — )■ g*, and let 
J'p = i*jR. For /i G g*, let be the stabilizer of p with respect to the coadjoint action, and consider 
the Marsden-Weinstein quotients {J^^{p)/Gp,ijj^)^ whose connected components are the symplectic 
leaves of tt^/gi see e.g. [2S]- Consider the natural embeddings 

I : JpHp)/G^. ^ Jr\p)/G,,. (A.81) 

Corollary A. 3. The symplectic leaves ofTTp/c are given by the connected components of J'f^ (p) / G ^ , 
for 1^1 Q* , with symplectic form l*uj^. 

Proof. The leaves of ttr/g are given by the connected components of J'p^{p)/Gf^, for ^ G g*, and by 
Lemma [A. H (ii), we have that the leaves of the connected components of {J'p ^(/u)/G^)n(P/G) 

and that they sit in {Jp^{pL) /G^) as symplectic submanifolds. Finally, note that 

{Jn\p)/G,) n (P/G) ^ iJpHo,) n P)/G = JpHo,)/g = Jp\p)/Gp, 

where is the coadjoint orbit through p. □ 

Suppose that -F is a regular distribution on P such that TP = F©Ker Op. The pair (F, 0,p) defines 
a bivector field vTp on P, as in (|2.10p . If F is G-invariant, then vTp is also G-invariant, and defines a 
bivector field on P/G via reduction. The relation between this bivector and the one in Lemma lA.l( i) 
is as follows. 

Lemma A. 4. The reduction of the almost Poisson structure Tip to P/G coincides with the Poisson 
structure ttp/q given by the reduction of the presymplectic structure Vtp as in Lemma \A.l\ i). 

Proof. We have to show that TTp^Q{a) = Tpp{'Kp{ppa)) for each 1-form a on P/G. Note that the vector 
field 7rl{ppa) satisfies i^j ^^,^^ilp|p = — ppa|p. So there is a 1-form F on P such that F|p = and 
it),, ,r2p = —ppO + T. Using that Ker^lp C Vp we see that i tt , * .^}p(X) = T(X) for all X G Ker rip 
and thus F = 0. Hence i^tt (^.^^^^p = —Pp<^, which implies that 7rp^Q(a) = Tpp{Trp{p*a)). □ 
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